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1 Introduction. 



New theories, new symmetries, new particles, new phenomena. 



The creation of the unified electroweak theory at the end of the 1960s |]], and of the 
quantum chromodynamics at the beginning of the 1970s || has dramatically changed the 
entire picture of elementary particle physics. Its foundation changed to gauge symmetries: 
electroweak symmetry SU(2)l x U(l) and the color symmetry SU(3) C . It became clear 
that the gauge symmetries determine the dynamics of the fundamental physical processes in 
which the key players are the gauge vector bosons — the well-familiar photon and a host of 
new particles: W + , W" , Z bosons and eight gluons which differ from one another in color 
charge. Even though the Higgs condensate, filling the entire space, remains enigmatic in the 
electroweak theory, while the problem of confinement is still unsolved in chromodynamics, 
the two theories are nevertheless so inseparable from modern physics that they were given 
the name of the Minimal Standard Model (MSM). We assume in this review that the reader 
is familiar with the basics of the MSM (see, for example, the monographs Q). 

The concept of quarks has undergone dramatic expansion in the process of creation of 
the MSM. In |l] the electroweak theory was suggested for leptons (electron and electron 
neutrino). The subsequent inclusion of quarks into the theory led to the hypothesis that in 
addition to the three quarks known at the time — u, d and s — there exists the fourth quark, 
c. According to 0, if d- and s-quarks are analogues, respectively, of e and fi, then the 
mutually orthogonal combinations ucos9 c + csin# c and —usin9 c + ccos9 c , where 9 C is the 
Cabibbo angle, must constitute the analogues of v e and ||. 

One of the important consequences of the electroweak theory was the prediction of the 
weak neutral currents. According to the theory, the neutral weak currents must be diagonal; 
in other words, neutral currents changing quark flavor (FCNC) are forbidden. This explained 
the absence of such decays as K° — > ee, K° — > ftfi and K + — > 7i + ee. Since there is no neutral 
current sd in the Lagrangian, these decays cannot occur in the tree approximation: they 
require loops with virtual H^-bosons. This is also true for the transitions K° — K° that are 
responsible for the mass difference between the K\ and i^-mesons. 

Diagonal neutral currents were discovered in reactions with neutrino beams ||, in rotation 
of the polarization plane of a laser beam in bismuth vapor @, and in the scattering of 
polarized electrons by deuterons ||. 

The charmed quark c was discovered in 1974 ||. Even before that, Kobayashi and 
Maskawa ]TD[ conjectured that in addition to two generations of leptons and quarks, (u e , e, u, d), 
(vn, fi, c, s), there must exist the third generation (u T , r, t, b). The 2x2 Cabibbo matrix for 
two generations, 

/ cos 0c sin 0c \ 
I — sin 9 C cos 9 C J 

is replaced in the case of three generations by a 3 x 3 unitary matrix that in its most 
general form contains three angles and one phase; the phase is nonzero if the CP-invariance 
is violated. This is how the mechanism of CP violation at the level of quark currents was 
proposed. 

The r-lepton |[l| and the 6-quark Jl2| were discovered experimentally in mid-1970s. The 
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Figure 1: The Z boson as a resonance in e + e~ annihilations. A fermion-antifermion pair in 
the Z-boson decay is denoted by //, where / can be either a lepton (e, fi, r, v e , v T ), or 
a quark (u, d, s, c, b). In this last case a quark-antiquark pair typically transforms, owing to 
interactions with gluons, to a multi-hadron state. The outgoing arrow in this and subsequent 
diagrams corresponds to the emission of the fermion (/) and to absorption of the antifermion 
(e + ); an incoming arrow denotes the emission of the antifermion (/) and the absorption of 
the fermion (e~). 



II PI. As for the 



heaviest fermion — the t-quark — was discovered only two decades later 
mechanism of CP-invariance violation, it remains unknown even now. 

The renormalizability of the electroweak theory and quantum chromodynamics ]15| and 
the property of asymptotic freedom in QCD [T(J opened a wide field for reliable computa- 
tions based on perturbation theory. On the basis of such computations in the tree-diagram 
approximation, it was possible to predict such qualitatively novel phenomena as quark and 
gluon jets; using the data on neutral currents, it proved possible to perform preliminary 
computations of the masses and partial widths of the W and Z bosons even before their 
actual discoveries. 



W and Z boson 'factories'. 



To test the predictions of the electroweak theory, proton-antiproton colliders were built 
at the beginning of the 1980s in Europe (at CERN) and then in the USA (at FNAL). The 
discovery of the W and Z bosons |L7| provided spectacular confirmation of the tree-diagram 
calculations M and made feasible and urgent precision tests of the electroweak theory with 
loops included. (When speaking about the tree (Born) approximation and loops, we always 
mean the corresponding Feynman diagrams.) 

A unique object for such tests was the Z boson. To carry out a precision study of its 
properties, electron-positron colliders were built at the end of the 1980s: LEP1 at CERN 
and SLC at SLAC. Electron and positron in these colliders collide at the center-of-mass 
energy equal to the Z boson mass. This results in the resonance creation and decay of 
a Z boson (see figure 1). The LEP1 completed its operations in October 1995; about 20 
million Z bosons were detected in the four detectors of this collider (ALEPH, DELPHI, L3, 
OPAL). The total number of Z bosons recorded at the SLC by the sole SLD detector was 
approximately 10 5 ; however, since the colliding electrons are longitudinally polarized, it was 
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Figure 2: Contribution of the tt to the Z-boson propagator (a). Contribution of the loop 
with virtual Z-boson and a higgs to the Z-boson propagator (b). 

possible to study the dependence of the annihilation cross section e + e~ of the Z boson on 
the sign of beam polarization. As a result, the SLC proved its competitiveness even with 
substantially lower statistics. The statistical and systematic accuracy achieved in the study 
of the properties of the Z-boson are of the order of 10 -5 for the Z boson mass and of the 
order of several thousandths for the observables that characterize its decays. 



What's the point of studying loop corrections? 



A natural question is: why do we need to compare the experimental data and the loop 
corrections of the electroweak theory? We need it mostly to gather data on the not yet 
discovered particles. For instance, even before the t-quark was discovered on the Tevatron 
by CDF |13j and DO [14| collaborations, its mass was predicted by analyzing the radiative 
loop corrections and the LEP1 and SLC data [IS|. The main loop involving the virtual t 



and t-quarks is shown in figure 2a. The scalar Higgs boson (or simply the higgs) had not 
been found yet. In the minimal version of the theory, the so-called Minimal Standard Model 
(MSM), there is a single higgs — a neutral particle whose mass is not fixed by the model. 
In the Minimal Supersymmetric Standard Model (MSSM) we have three neutral and two 
charged higgses. The lightest of the neutral higgses must not be heavier than 135 GeV [ |T9|1 , 



2f| . The simplest diagram involving a virtual higgs is shown in figure 2b. 

When planning experiments on LEP1 and SLC, people had great expectations that preci- 
sion measurements would detect pronounced deviations from the predictions of the standard 
model and would thereby unambiguously point to the reality of some sort of 'new physics'. 
In fact, even though some discrepancies with the MSM were found, they go beyond the 
three standard deviations only in a single case (that of the decay of a Z boson to a 66-pair). 
If these discrepancies are not caused by some sort of systematic error, they may indicate 
(see Conclusions) the existence of the relatively light (~ 100 GeV) squarks and gluino — the 
supersymmetric partners of quarks and gluons, respectively. 
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2 Brief history of electroweak radiative corrections. 



The pioneer calculations of electroweak corrections in MSM were performed in the 1970s, 
long before the discovery of the W and Z bosons. The calculations were devoted to the 
muon decay and the /3-decays of the neutron and nuclei and to deep inelastic processes. In 
connection with the construction of LEP and SLC, a number of teams of theorists carried out 
detailed calculations of the required radiative corrections. These calculations were discussed 
and compared at special workshops and meetings. The result of this work was the publication 



of two so-called 'CERN yellow reports' [^TJ, | 22| , which, together with the 'yellow report' 



| T8|j , became the 'must' books for experimenters and theoreticians studying the Z boson. 

Muon and neutron decays. 



Sirlin |23| calculated the radiative corrections to the muon decay due to one-loop Feynman 
diagrams (see figure 3). We must emphasize that the purely electromagnetic correction to 
the muon decay due to the exchange of virtual photons and the emission of real photons was 
calculated even earlier [24j| , for the pointlike four-fermion interaction, i.e. without taking the 
H^-boson into account (see figure 4). It was found that the correction is finite: it contained 
no divergences. The four-fermion interaction constant G M , extracted from the muon lifetime 



" 192tt 3 77?.^ 2vr 4 V ' 

where f(x) = 1 — 8x + 8x 3 — x 4 — 12x 2 lnx, already includes this electromagnetic correction 
proportional to a. G M = 1.16639 (2) \0~ 5 GeV~ 2 . The finiteness of the purely electromagnetic 
correction in the muon decay is caused by the V-A nature of the interacting charged currents 
u^fi and ev e . 

In the neutron decay, the purely electromagnetic correction to the four-fermion interac- 
tion (figure 4(j)) diverges logarithmically. In view of the W^-boson propagator (figure 4(h)), 
the logarithmic divergence is cut off at the H^-boson mass. 

This correction to the vector vertex in the leading logarithmic approximation, calculated 

in 



25], is given by the factor 



1 + ^ln^, (2) 

27r m p 

where m p is the proton mass. Numerically, its value is of the order of 1.7%. Only after 
this correction is taken into account, does cos 2 9 C extracted from the nuclear /3-decay become 
equal to 1 — sin 2 9 C , where sin 2 8 C is found from the decays of strange particles (9 C is the 
Cabibbo angle). Although the correction we discuss now contains mw in the logarithmic 
term, it is essentially electromagnetic and not electroweak, since it is insensitive to details 
of the electroweak theory at short distances, in contrast to, say, electroweak corrections to 
the muon decay (figure 3). 

Calculations of electroweak corrections to the muon decay show that the main contribu- 
tion, exceeding all others, is caused by the vacuum polarization of the photon (figure 5a). 
At the first glance, this correction should not have emerged in the /z-decay in the one-loop 
approximation: it is not there among the loops in figure 3. It does appear, however, when 
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Figure 3: One- loop diagrams in muon decay. Loops in the VF-boson propagator (a), (b), 
(c), (d); in the W^-vertex (e), (f) and (g) and in external fermion lines (h), (i), (j), (k) (similar 
diagrams for e and V e are assumed), as well as square- type diagrams (1), (m), (n), (o). 
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a) b) 




Figure 4: Muon decay in tree-diagram approximation (a) and in local four-fermion approx- 
imation (b). Electromagnetic corrections to muon decay in the local approximation (c) -(g). 
Electromagnetic corrections to the /3-decay of the neutron: in tree-diagram approximation 
(h), in four-fermion approximation (j). 
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Figure 5: Photon polarization of the vacuum, resulting in the logarithmic running of the 
electromagnetic charge e and the 'fine structure constant' function of q 2 , where 

q is the 4-momentum of the photon (a). Some of the diagram that contribute to the proper 
energy of the W^-boson (b)-(g). Some diagrams that contribute to the proper energy of the 
Z-boson (h)-(n). Some of the diagram that contribute to the Z <-> 7 transition (o)-(r). 
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G M is expressed in terms of the fine structure constant a and the masses of the W and Z 
bosons. 



Main relations of electroweak theory. 



It is well known ||] that in the Born approximation we have (see figure 4a) 



where mw is the W boson mass and g is its coupling constant to the charged current. 
On the other hand, we have in the same approximation 

m w = ^9V , ( 4 ) 

where 77 is the vacuum expectation value of the higgs field. 
Likewise, 

m z = ^fv , ( 5 ) 

where mz is the Z boson mass and / is its coupling constant to the neutral left-handed 
current. 

Therefore, 

m w 9 / fi x 

= 7 • (6) 

m z } 

If we introduce the famous Weinberg angle flT], it becomes obvious that the two definitions 
are equally valid in the Born approximation: 

m w 

costV = , 

m z 

cos6 w = j. (7) 

It is well known |l|] that the angle 9w in the electroweak theory defines the relation between 
the electric charge e and the weak charge g: 

e = g sin 9 W . (8) 

In the Born approximation, therefore, 

G g2 1 — = — . (9) 

M 4A/2m^ V^m^sin 2 ^ y/2m 2 z sin 2 9 W cos 2 9 W 

If we now take into account the electroweak corrections, we find 

m w ^ g_ 
m z f 
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Traditional parametrization of corrections to the /i-decay and the running a. 



Sirlin's definitions |23[] are widely used in the literature (see the review by Langacker and 
Erler [^] and references therein); according to them 

2 

= sin 2 6 W = 1 - c 2 ^ = 1 - cos 2 9 W = 1 \ , (10) 

m z 

na 

" v^mf^(l-Ar) ' 1 J 

where 

Ar = Ar em + Ar eu , 

includes both the truly electroweak correction Ar ew for loops in figure 3, and the purely 
electromagnetic correction Ar em due to a running from q 2 = to q 2 ~ myy,m%. This 
correction arises because 

a = a(q 2 = 0) = 1/137.035985(61) (12) 

is defined for q 2 = 0, while typical momenta of virtual particles in the electroweak loop are 
of the order of the intermediate boson masses. It is convenient to denote 

a = a(m 2 z ) = ^ = " . (13) 

1 - Ar Rm 1 - da 



Ar em was calculated in a number of papers |2(J ; the necessary formulas are given in Appen- 
dices B and C. The contribution of the lepton loops to Ar em is described by the expression 

5a l = Ar l em = f J2[\n ™§ - |] = 0.03141 , (14) 

where I — e,fi, t. The contribution of quark loops cannot be calculated theoretically because 
the quark mass in the logarithm has no rigorous theoretical definition. This reflects our 
ignorance of strong interactions at large distances (small momenta). The quark (hadron) 
part Ar-g m is therefore calculated by substituting into the dispersion relation the experimental 
data on the cross section of the e + e~-annihilation into hadrons: 

6a h = Ar h em = ^f- f°° , (15) 

where cjg +e _ is the cross section of the e + e~ annihilation into hadrons via one virtual photon. 

In this review we make use of the recent result reported by Jegerlehner and Eidelman 
2§: Ar h em = 0.02799(66), so that 



5a = 5a l + 5a h = Ar em = Ar l em + Ar h em = 0.05940(66) . (16) 
As follows from ( |i~3"|) and ([16]), 



a 



[128.896(90)]" 1 . (17) 
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A summary of results of various calculations of a is given in Appendix C. Following tradition, 
the contributions of the t-quark loop and the W-boson loop are not included into a(m^). 
In the leading approximation in l/mf, the contribution of the t-quark loop is 

Ar* = f = -0.00005 for m t = 180GeV , (18) 

6771 7r 45 m t 

and the exact formula (see equation (^)) corresponds to 

Arl m = -0.00006 for m t = 180GeV . 

The contribution of the W^-boson loop is gauge- dependent. In the 't Hooft-Feynman gauge 
it is Ar e ^ = 0.00050 (see <M)). 



Deep inelastic neutrino scattering by nucleons. 

A predominant part of theoretical work on electroweak corrections prior to the discovery 
of the W- and Z-bosons was devoted to calculating the neutrino-electron [2~7| and especially 



nucleon-electron [28| interaction cross sections. The reason for this is that after the discovery 



of neutral currents the quantity = 1 — -^f 1 was extracted precisely from a comparison of 

the cross section of neutral currents (NC) and the charged currents (CC). While a M^-boson 
interacts with a charged current of V — A type, for example, 

^W a u{ la + lal5 )d , (19) 

the Z-boson interacts with neutral currents that have a more complex form, 

f -Z a ^ f [Ti lal5 + {T( - 2Qh 2 w ) la }^ f , (20) 

where T% is the third projection of the weak isotopic spin of the left-hand component of 
the fermion / (quark or lepton), Q* is its charge, and ipf is the Dirac spinor describing it. 
In the tree approximation, the NC/CC cross section ratio for purely axial interactions (and 
isoscalar target) equals 1, since in this approximation the quantity 

P = 4 ! 4 (21) 

equals unity. With the vector current taken into account, the ratio of NC and CC is a 
function of s\y- Measurements of this ratio gave w 0.23, which thus made it possible to 
predict the masses of the W and Z bosons (from the formula for the muon decay @). More 
accurate measurements of NC/CC made it possible to improve the accuracy of to such 
an extent that it became necessary to take into account the electroweak radiative corrections 
both in Syy and in p, which now, with the corrections taken into account, is not equal to 
unity any more. 

Veltman was the first to point out [^jj that if m t /mz 1, the main correction to p 



and s^r is caused by the violation of the electroweak isotopic invariance by the masses of 
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the t- and fe-quarks in loops of self-energies of the Z and W^-bosons. To find p in the limit 
m t /mz 1, it is sufficient to consider these loops neglecting the momentum of W and 
Z bosons q in comparison with the masses of the W and Z, i.e. for q 2 = 0. Elementary 
calculation of the loops indicated above yields (see Appendix H): 

p=j | 3az ( mt )2 _ 1 | 3a ^( m t Y - I | 3G ^ m i (22) 



167T m z 167T m w 8\/2' 



7T 



Here and hereafter we denote 



P g 2 

olz = — , ol w = — . (23) 

47T 47T 

Since in real life m t /mz — 2, the sum of the remaining, non-leading corrections is found 
to be comparable to the correction proportional to m\ (see below). 

At the present moment, the quantity Sw extracted from the data on deep inelastic vN- 
scattering is determined as 0.2260(48) (the global fit of the data from the collaborations 
CDHS [0 and CHARM |3T|), or 0.2218(59) (collaboration CCFR) [0. The accuracy 



of these data is poorer than found by the direct measurement of the W boson mass mw 
and of the ratio mw/mz by the collaborations UA2 in CERN []33j and CDF at Tevatron 
34|. According to PDG |55j, the fitted quantity is mw = 80.22(26), which corresponds 



to s w = 0.2264(25). Note that according to the most recent data p5| , the measurement 
accuracy is even higher for mw- — 80.26(16)Gel / . For this reason we will not discuss 
further the deep inelastic scattering of the neutrino. Moreover, additional assumptions on 
the effective mass of the c-quark and on the accuracy of the QCD corrections are necessary 
for the interpretation of these experiments. 

Other processes involving neutral currents. 

This is true to even greater degree for the parity violation in eD-scattering ||, which, in 
addition, provides a considerably less accurate s w = 0.216(17). 

The measurement accuracy of the u^e- and z/^e-scattering even in the highest-accuracy 



experiment CHARM II |37] is not sufficient for revealing the genuine electroweak corrections. 



However, after an analysis in the Born approximation, this experiment has demonstrated 
for the first time that the interaction constant of the current v^v^ with the Z boson is in 



satisfactory agreement with the theory (see |38| , [35] ) . 



The experiment on measuring parity violation in cesium 133 Cs55 [^0f is also insufficiently 
sensitive, at the accuracy achieved. Here the effect is produced by the interaction of the 
nucleon vector current with the electron axial current. Since the characteristic momentum of 
electron in an atom is small in comparison with nuclear dimensions, all nucleons of a nucleus 
'function' coherently and the nucleus is characterized by an aggregate weak charge Qw, 
which is experimentally found to be Q W P — —71.0 ± 1.8, while the theoretically anticipated 
value is —72.9 ± 0.1. A spectacular property of Qw is the fact that owing to an accidental 
cancellation of protons' and neutrons' contributions, it is practically independent of m t . On 
the contrary, Qw is very sensitive to the contribution of neutral bosons Z' and Z" heavier 
than the Z-bosons (if they exist). 

The best object for testing the electroweak theory at the loop level is therefore the Z 
boson with which this review is predominantly connected. 
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3 On optimal parametrization of the theory and the 
choice of the Born approximation. 



The electroweak theory is in many ways different from electrodynamics, in particular in the 
diversity of particle interactions that must be taken into account when considering any effect 
in the loop approximation. The parametrization of QED is straightforward: fundamental 
quantities are the electron mass and charge, which are known with very high accuracy. It is 
therefore natural to express all theoretical predictions of QED in terms of a and m e . 

Traditional choice of the main parameters. 

The choice of the main parameters in electroweak theory is not equally obvious. His- 
torically, those selected were as an experimentally measured with best accuracy weak 
decay constant, s 2 ^ = 1 — m^/m|, since W and Z bosons were not yet directly observed 
at that time, while the value of sfy was known from experiments with neutral currents, and 
finally, a. This parametrization of the 1970s proved to be surprisingly long-lived; the loop 
parameters Ar and p connected with it are widely used in the literature and are very likely 
to survive beyond the end of this century. 

In fact, this parametrization is far from being optimal because (and thus syy as well) 
is measured experimentally at much poorer accuracy than mz' 

m w = 80.26(16)GeV (24) 

m z = 91.1884(22)GeV (25) 

As a result, s\y is extracted by fitting the loop formulas for various observables. This 
extraction inevitably requires that we fix the values of the t-quark mass and the higgs 
mass. Another drawback of this parametrization is that the quantity a, despite its superior 
accuracy, is not directly related to the electroweak loops, which are characterized by a 
quantity a, and these we know with much less impressive accuracy. As a result, the purely 
electromagnetic correction Ar em is not separated from the genuinely electroweak corrections, 
thus blurring the interpretation of the experimental data. 

Optimal choice of the main parameters. 



As follows from the above remarks, the currently optimal parametrization is the one 
based on G M , m z and a. With this parametrization it is convenient to introduce the weak 
angle 9, defined (by analogy to equation (§)) by the relation 



where s 2 = sin 2 9, c 2 = cos 2 9 
As follows from equation (E 



M ^/2m 2 z s 2 c 2 ' ^ ^ 



sin 2 29 = —= = 0.71078(50) 



V2G, 
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s 2 = 0.23110(23) , (27) 
c = 0.87687(13) . 

The angle 9 was introduced in mid-1980s However, its consistent use began only 
after the publication of [[Kj. 

Using 9 instead of 9w automatically takes into account the running of a and makes 
it possible to concentrate on the genuinely electroweak corrections. Using mz instead of 
sw allows one to explicitly single out the dependence on m t and ran f° r eacri electroweak 
observable. 

Note that a different definition of the Z boson mass m z is known in the literature, related 
to a different parametrization of the shape of the Z boson peak J43|. This mass mz is smaller 



than mz by approximately 30 MeV. In this review, we consistently use only mz, following 



the summary reports of LEP collaborations [44]. 



Let us show how the parametrization in terms of G^, mz and s is applied to the decay 
amplitudes of Z boson and to the ratio m w /m z . 

Z boson decays. Amplitudes and widths. 



In correspondence with equation (p0j) , we rewrite the amplitude of the Z boson decay 
into a fermion-antifermion pair // in the form 

M(Z - ff) = ^fZ a $ f { la g Vf + lo nsg A Mf • (28) 

Here by definition / is the value of the coupling constant / in the Born approximation, 

j 2 = 4 v / 2G>i| = 0.54866(4) (29) 

(the use of the same letter / to denote both the fermion and the coupling constant cannot 
lead to confusion, to such an extent are these objects different). The high accuracy with 
which the numerical value of / is known comes from the fact that / is independent of 
a. All electroweak radiative corrections are 'hidden' in dimensionless constants gvf and 
QAf- These coefficients do not include the contribution of the final state interactions due to 
the exchange of gluons (for quarks) and photons (for quarks and leptons). The final state 
interactions have nothing in common with the electroweak corrections and must be taken 
into account as separate factors in the expressions for the decay rates. These factors are 
sometimes known as 'radiators', since they cover not only exchange of photons and gluons 
but also their emission. 

Radiators are trivially unities in the case of decay to any of the neutrino pairs v e v e , v^v^, 
v T v T and therefore 

T v = T{Z^vV)=±T Q {g Av + g 2 Vv ) , (30) 
where Tq is the so-called standard width: 

G 777,^ 

T = "/ = 82.944(6)MeV. (31) 
24y 2tt 
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If neutrino masses are assumed to be negligible, then 



9 Ay = 9v„ = 9u (32) 

so that 

r„ = ST gl. (33) 
For decays to any of the pairs of charged leptons 11 we have [18 



r(z - ii) = 4r [^(i + 2°) + 9 l(i + ^ - e4)l , (34) 

47T 47T mc 



where the QED correction is taken into account only to the lowest approximation in a; we 
neglect terms on the order of (a/71) 2 ~ 10~ 6 . The term proportional to mf is negligible for 
I — e,fi and must be included only for I = r {w^./m 2 z = 3.8 ■ 10~ 4 ). 
For decays to any of the five pairs of quarks qq we have 

T q = T(Z - qq) = 12 • T [g 2 Aq R Aq + g 2 Vq R Vq ]. (35) 

Here the factor 3, additional in comparison with leptons, takes into account the three colors 
of each quark. The 'radiators' in the first approximation are identical for the vector and the 
axial- vector interactions, 

a 

R Vq = R Aq = 1 + , (36) 

71 

where a s is the constant of interaction of gluons with quarks at q 2 = m|. There are different 
conventions for the choice of a s . In calculations of decays of the Z-boson, it is quite typical 
to determine a s using the so-called modified minimal subtraction scheme, MS (see at the 
end of Appendix A). The numerical value of a s , found from Z-boson decays, is of the order 
of 0.12. For additional details on the value of a s and for more accurate expressions for 
radiators see Appendix F. Here we only remark that vast literature is devoted to radiator 
calculations; they are calculated using perturbation theory up to terms (a s /7i) 3 , see [fR|, 



46], [47], [48]. The full hadron width is, to the accuracy of very small corrections, the sum 



of widths of five quark channels: 

r h = r u + r d + r s + r c + r fe . (37) 

The full width of the Z boson is given by the obvious expression: 

r z = r h + r e + r /1 + r T + 3r, . (38) 

The annihilation cross section e + e~ into hadrons at the Z-peak is given by the Breit-Wigner 
formula: 

i27rr e r, 

* h = mTf ■ (39) 

Finally, the following notation for the ratio of partial widths is widely used: 

R b = ^ , R c = ^ , Ri = ^ ■ (40) 

1 h 1 h t I 

(Note that in contrast to Rb and R c , Th in Ri is in the numerator.) 
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Asymmetries. 



In addition to the total and partial widths of Z boson decays, experimentalists also 
measure effects due to parity non-conservation, i.e. the interference of the vector and axial- 
vector currents. For pairs of light quarks (u, d, s, c) and leptons we determine the quantity 

Af = \ 9Ai9V l ■ (41) 
9Af + 9v f 

For the pair bb (see Appendix G), 

A 2gAb9vb (A9 , 

b ~ v*g Ab + i(3-v*) g z vb - V ' {A2} 

where v is the velocity of the 6-quark (in units of c): 



\ 



(43) 



Here rhb is the value of the 'running mass' of the 6-quark with momentum m^, calculated in 



the scheme MS ||50|| . The forward-backward charge asymmetry in the decay to // equals 
(see Appendix G) 

. t Np — Nb 3 

s N^N~ B = I A ' A > ■ (44) 

where A e refers to the creation of a Z boson in e + e~-annihilations, and Af refers to its decay 
into //. 

The longitudinal polarization of the r-lepton in the Z — > rf decay is 

P T = A T . (45) 

If, however, we measure polarization as a function of the angle 9 between the momentum 
of a t~ and the direction of the electron beam, this permits the determination of not only 
A T but A e as well: 

A T (l + cos 2 fl) + A e 2cosfl 
PAC0S6) = - (l + cos'g) + ^ e 2cosg • (46) 

The polarization of P T is found from P T (cos6) by separately integrating the numerator 
and the denominator in (|4q) over the total solid angle. 

The relative difference between total cross sections in the Z peak for the left- and right- 
polarized electron that collide with non-polarized positrons (this quantity is measured at the 
SLC collider) is 

Alr ss ^5 = A e . (47) 

The measurement of the asymmetries outlined above allows one to experimentally de- 
termine the quantities gvf/dAf, since these asymmetries are caused by the interference of 
the vector and axial-vector currents. In their turn, measurements of the widths Tf, T^, Tz 
mostly permit the experimental determination of gAf, since \gvq\ 2 < \gA q \ 2 for quarks, and 
for leptons \gvi\ 2 <^ |fiU/| 2 - As for T q and Th, getting these quantities allows one to find a s . 
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The Born approximation for hadronless observables. 



Before discussing the loop electroweak corrections, let us consider expressions for mw/mz, 
gAf and gvf/dAf i n the so-called a-Born approximation. Using the angle 9 introduced earlier, 
its sin 9 = s and cos 9 = c, we automatically take into account the purely electromagnetic 
correction due to the running of a. It is easily shown that in the a-Born approximation 

(m w /m z ) B = c (48) 

9A f = T 3f (49) 
(g Vf /9Af) B = i-MQf\s 2 (50) 

It is of interest to compare the Born values with their experimental values. Table 1 presents 
this comparison for the so-called 'hadronless' observables (here and below the experimental 
results are taken from f44j.) 



Table 1 



Observable 


Experiment 


o-Born 


m w /m z 


0.8802(18) 


0.8769(1) 


m w (GeV) 


80.26(16) 


79.96(2) 


b w 


0.2253(31) 


0.2311(2) 


9AI 


-0.5011(4) 


-0.5000(0) 


Ti (MeV) 


83.93(14) 


83.57(2) 


Qvi/qai 


0.0756(14) 


0.0756(9) 


s 2 


0.2311(4) 


0.2311(2) 



For the reader's convenience, the table lists different representations of the same observables 
known in the literature. Thus, according to widely used definitions, 

s 2 w = 1 - m 2 w /m 2 z , (51) 

sf = s 2 eff = sin 2 <%} = 1(1 - g vl /g Al ) . (52) 

The experimental value of sf in the table is the average of two numbers, 0.2316(5) (LEP) 
and 0.2305(5) (SLC). It is assumed in the table that the lepton universality holds, thus the 
lepton decay data have been averaged over a number of observables. 

Table 1 shows that the a-Born approximation provides good description of the exper- 
imental data. ( An agreement is found for the hadron decays of Z-bosons as well). An 
especially (and unexpectedly!) good agreement, and the ensuing smallness of radiative cor- 
rections, is found for gvi/gAi- The anomalous smallness of true electroweak corrections was 
first pointed out in 1992 f4~2| , when the a-Born approximation was applied for the first time. 
(Before that it was hidden in the shadow of the large contribution of purely electromagnetic 
running of a, that was not separated from truly electroweak corrections). The LEP1 data 
of 1992 were not sufficiently accurate to allow detecting them. Even the data presented to 
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the Marseille conference in summer 1993 were not, as pointed out in a number of reports 
5T], sufficiently accurate for the detection of radiative corrections. At the Glasgow confer- 



ence in summer 1994 corrections were detectable at the level of 2.3a for g^i, 1.5a for mw 
and la for gvi/gAi Note that the difference between most recent experimental values 



s w = 0.2253(31) and sf = 0.2311(4) is a 2a manifestation of electroweak radiative correc- 
tions, which is independent of either the choice of the Born approximation or of the choice 
of calculation scheme. 



4 One-loop corrections to hadronless observables. 

The fact that the experimentally observed electroweak corrections are small suggests that 
one-loop approximation would be sufficient for describing them [ 4"2"| . This idea is supported 



by a thorough evaluation |L8[ of theoretical uncertainties contributed by higher-order pertur- 
bation theory in electroweak interaction. There exists essentially a single two-loop diagram 
that should be taken into account. It was calculated in |f>3|l, and we will discuss its contri- 



bution and take it into account below (see formula ([79])). 

Four types of Feynman diagrams. 



Four types of Feynman diagrams contribute to electroweak corrections for the observables 
of interest to us here, m w /m z , 9ai, 9vi/9ai- 

1. Self-energy loops for W and Z bosons with virtual z/, /, q, H, W and Z in loops. Exam- 
ples of some of these diagrams are shown in figures 5b-5n. 

2. Loops of charged particles that result in Z <-> 7 transitions (figures 5o-5r). 

3. Vertex triangles with virtual leptons and a W or Z boson (figures 6a-6c). 

4. Electroweak corrections to lepton wavefunctions (figures 6d-6e). 

It must be emphasized that loops shown in figures 5h-5n contribute not only to the mz 
mass and, consequently, to the mw/mz ratio but also to the Z boson decay to U, to which 
Z <-> 7 transitions also contribute (figures 5o-5r). This occurs owing to the diagrams of the 
type of figures 6f-6g which give corrections to the Z-boson wavefunction. 

Obviously, electroweak corrections to m w /m z , Qai and gvi/dAi are dimensionless and 
thus can be expressed in terms of a, c, s and the dimensionless parameters 

t = (^) 2 , /> = (— ) 2 , 
m z m z 

where m t is the mass of the t-quark and m# is the higgs mass. (We neglect the masses of 
leptons and all quarks except t.) 

The asymptotic limit at m^ 3> m|. 
The functions V m (t,h), V A (t,h) and V R (t,h). 
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1 1 



a) 



b) 



c) 





Figure 6: Vertex triangular diagrams in the Z — > 11 decay (a), (b), (c). Loops that renormal- 
ize the lepton wavefunctions in the Z — > 11 decay. (Of course, antilepton have similar loops.) 
(d), (e). Types of diagrams that renormalize the Z boson wavefunction in the Z — > 11 decay 
(f), (g). Virtual particles in the loops are those that we were discussing above 
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It is convenient to split the calculation of corrections into a number of stages and begin 
with calculating the asymptotic limit for t ^> 1. 

According to the reasons mentioned above (see equation fl2~2"|) , the main contribution 
comes from diagrams that contain t- and 6-quarks (5c and 5i,j). A simple calculation (see 
Appendix H) gives the following result for the sum of the Born and loop terms: 

3c 

mw/mz = C+ 32*s*(<?-s*) m ■ (53) 
1 3 

9Ai = --~ — Y^ at ' 54 

3 

R = gvi/gAi = 1 - 4s 2 + 4?r(c2 _ s2) «t , (55) 

1 3 

2 b4irs z c z 

The functions V m (t,h), V A (t,h) and V R (t,h). 



If we now switch from the asymptotic case of t ^> 1 to the case of t ~ 1, then, first, the 
change in the contribution of the diagrams 5c, 5i and 5j can be written in the form 

t->* + T,(t) , (57) 

where the index i = m, A, R, v for m w /m z , gAi, R = Qvil 9ai and g u , respectively. 

The functions Tj are relatively simple combinations of algebraic and logarithmic functions. 
They are listed in explicit form in Appendix I. Their numerical values for a range of values 
of mt are given in Table 2. The functions Tj(t) thus describe the contribution of the quark 
doublet t, b to mw/mz, gA, R = gvi/gAi and g v . If, however, we now take into account the 
contributions of the remaining virtual particles, then the result can be given in the form 

t^V i (t,h) = t + T i (t) + H i (h) + C i + 6Vi(t,h) . (58) 
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Table 2 



m t 
(GeV) 


t 


T 


T A 


Tr 








-0 


188 


0.875 


0.444 


10 


0.012 





192 


0.934 


0.038 


20 


0.048 


-0 


256 


0.955 


-0.015 


30 


0.108 


-0 


430 


0.812 


-0.305 


40 


0.192 


-0 


753 


0.403 


-0.959 


50 


0.301 


-0 


985 


0.111 


-0.748 


60 


0.433 


-0 


931 


0.327 


-0.412 


70 


0.589 


-0 


688 


0.390 


-0.250 


80 


0.770 


-0 


317 


0.421 


-0.143 


90 


0.974 


-0 


080 


0.440 


-0.061 


100 


1.203 





084 


0.451 


0.006 


110 


1.455 





214 


0.460 


0.062 


120 


1.732 





323 


0.465 


0.111 


130 


2.032 





418 


0.470 


0.154 


140 


2.357 





503 


0.473 


0.193 


150 


2.706 





579 


0.476 


0.228 


160 


3.079 





649 


0.478 


0.261 


170 


3.476 





713 


0.4^ 


30 


0.291 


180 


3.896 





772 


0.4^ 


31 


0.319 


190 


4.341 





828 


0.4^ 


33 


0.345 


200 


4.810 





880 


0.4^ 


34 


0.370 


210 


5.303 





929 


0.4^ 


35 


0.393 


220 


5.821 





975 


0.4^ 


35 


0.415 


230 


6.362 


1 


019 


0.4^ 


36 


0.436 


240 


6.927 


1 


061 


0.4^ 


37 


0.456 


250 


7.516 


1 


101 


0.4^ 


37 


0.475 


260 


8.130 


1 


139 


0.4^ 


37 


0.493 


270 


8.767 


1 


176 


0.4^ 


38 


0.510 


280 


9.428 


1 


211 


0.4^ 


38 


0.527 


290 


10.114 


1 


245 


0.4^ 


39 


0.543 


300 


10.823 


1 


277 


0.4^ 


39 


0.559 



Here Hi(h) contain the contribution of the virtual vector and higgs bosons W, Z and H 
and are functions of the higgs mass m# . (The masses of the W and Z bosons enter Hi(h) 
via the parameters c, s, defined by equation (p6|)). The explicit form of the functions Hi is 
given in Appendix I, and their numerical values for various values of ran are given in Table 
3. 
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Table 3 



m H 
(GeV) 


h 


H m 


H A 


Hr 


0.01 





000 


1.120 


-8.716 


1.359 


0.10 





000 


1.119 


-5.654 


1.354 


1.00 





000 


1.103 


-2.652 


1.315 


10.00 





012 


0.980 


-0.133 


1.016 


50.00 





301 


0.661 


0.645 


0.360 


100.00 


1 


203 


0.433 


0.653 


-0.022 


150.00 


2 


706 


0.275 


0.588 


-0.258 


200.00 


4 


810 


0.151 


0.518 


-0.430 


250.00 


7 


516 


0.050 


0.452 


-0.566 


300.00 


10 


823 


-0.037 


0.392 


-0.679 


350.00 


14 


732 


-0.112 


0.338 


-0.776 


400.00 


19 


241 


-0.178 


0.289 


-0.860 


450.00 


24 


352 


-0.238 


0.244 


-0.936 


500.00 


30 


065 


-0.292 


0.202 


-1.004 


550.00 


36 


378 


-0.341 


0.164 


-1.065 


600.00 


43 


293 


-0.387 


0.128 


-1.122 


650.00 


50 


809 


-0.429 


0.095 


-1.175 


700.00 


58 


927 


-0.469 


0.064 


-1.223 


750.00 


67 


646 


-0.506 


0.035 


-1.269 


800.00 


76 


966 


-0.540 


0.007 


-1.311 


850.00 


86 


887 


-0.573 


-0.019 


-1.352 


900.00 


97 


410 


-0.604 


-0.044 


-1.390 


950.00 


108 


534 


-0.633 


-0.067 


-1.426 


1000.00 


120 


259 


-0.661 


-0.090 


-1.460 



The constants Ci in ( |58D include the contributions of light fermions to the self-energy of the 
W and Z bosons, and also to the diagrams of figure 3, describing the muon decay, and of 
figures 6a-6c, describing the Z-boson decay. The constants Ci are relatively complicated 
functions of s 2 (see Appendix M). We list here their numerical values for s 2 = 0.23110 — 5s 2 : 



C m = -1.3497 + 4.135s 2 , 


(59) 


C A = -2.2621 - 2.635s 2 , 


(60) 


C R = -3.5045 - 5.725s 2 , 


(61) 


C u = -1.1641 -4.885s 2 . 


(62) 


Corrections 5V ; (t,h). 





Finally, the last term in equation (|58|) includes the sum of corrections of five different 
types. Their common feature is that they are all quite small (except for 5\Vj) and that they 
represent two loops (with the exception of a one-loop 8{Vi and a three-loop 8$Vj). 
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t 

a) 




c) 



Figure 7: Virtual t-quarks (a) and W bosons (b), (c) in the photon polarization of the 
vacuum. 



5\Vi contains contributions of the PF-boson and the t-quark to the polarization of the 
electromagnetic vacuum 8 w a = Ar^ and 8 t a = Ar* m (see figures 7a-c and equation 
(|16D), which traditionally are not included into the running of ct(q 2 ), i.e. into a. It 
is reasonable to treat them as electroweak corrections. This is especially true for the 
iy-loop that depends on the gauge chosen for the description of the W and Z bosons. 
Only after this loop is taken into account, the resultant electroweak corrections become 
gauge-invariant, as it should indeed be for physical observables. Here and hereafter in 
the calculations the 't Hooft-Feynman gauge is used (see Appendix A), 



where 



16 1 

6iV m (t, h) = 7ts 4 — {5wct + 5 t a) = —0.055 

3 a 

16 1 

SxVrU, h) = ns 2 c 2 -(5 w a + S t a) = -0.181 

3 a 

5 1 V A {tih)=5 1 V u {t,h) = U, 



J*^ = -L[(3 + 4c 2 )(l - V4c 2 - 1 arcsin ^-) - \\ = 0.0686 , 
a Z7T 2c 3 

8 t a 4 r . , . , l n 4 1 
~~ = -Q- 1 + 2t)F t (t) - - ~ - — - + ... ~ -0.00768 . 
a 97r 3 457T t 



(63) 

(64) 
(65) 

(66) 
(67) 
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(Here and in what follows, unless specified otherwise, we use m t 
ical evaluations.) 



175 GeV in numer- 



2. The corrections b^Vi are caused by including virtual gluons in electroweak loops in the 
order aa s (see figures 8a-c); similar diagrams can, of course, be drawn for H^-bosons. 
In addition to loops with light quarks q = u,d, s, c, there exist similar loops with 
third-generation quarks t and b: 

5 2 V i (t) = 8$V i + 5 t 2 V i (t) . 

The analytical expressions for corrections S^Vi and 8^Vi(t) are given in Appendix O. 
Here we only give numerical estimates for them, 



fi q V 



-0.377- 



71 



5lV A 



1.750- 



.a. 



71 



51V R = 



slv m {t) 



-11.67 



5lV A (t) = -10.10 
%V R (t) = -11.88 



<y s {m t ) 

71 

a s (m t ) 



7Y 

a s (m t ) 



-10.61- 

7T 

do 

-9.18— 

7T 



-10.80- 



7Y 



7Y 



since 



16 1 



1 + j§-a s hit 

12tt " 



(68) 

(69) 
(70) 

(71) 
(72) 
(73) 

(74) 



(For numerical evaluation, we use a s = a s (mz) = 0.125.) We have mentioned already 
that the corrections S^V^t), whose numerical values were given in Q7TD— (J73[) , are much 
larger than all other terms included in SVi. We emphasize that the term in 5 l 2 Vi that 
is leading for high t is universal: it is independent of i. As shown in [53], this leading 
term is obtained by multiplying the Veltman asymptotics t by a factor 

2tt 2 + 6 a s (m t ) 



1 - 



9 



7T 



or, numerically, 



t -+t(l -2.86 



a s (m t ). 
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(75) 



(76) 



Qualitatively the factor (ffq) corresponds to the fact that the running mass of the t- 
quark at momenta p 2 ~ mf that circulate in the t-quark loop is lower than "on the 
mass-shell" mass of the t-quark. It is interesting to compare the correction ( [76]) with 
the quantity 



2/ 



-m. 



2.78 



a s (m t ). 



71 



(77) 



calculated in the Landau gauge in ||55|| , p 102. The agreement is overwhelming. There 
is, therefore, a simple mnemonic rule for evaluating the main gluon corrections for the 
i-loop. 
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a) * b) 




c) d) 




H 



j) 

Figure 8: Gluon corrections to the electroweak quark loop of the proper energy of the Z 
boson (a)-(c). Higgs corrections to the electroweak t quark loop of the Z boson (d)-(f). 
Two-loop higgs corrections (g)-(j). 
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3. Corrections 5^Vi of the order of aa 2 were calculated in the literature |56| for the term 
leading in t (i.e. aa 2 t). They are independent of i (in numerical estimates we use for 
the number of quark flavors Nf = 5): 

5 3 Vi(t) ~ -(2.38 - 0.l8N f )a 2 (m t )t ~ -lA8a 2 (m t )t = -0.07 . (78) 

The corrections SiVi, 5 2 Vi 1 S 3 Vi, are independent of m#, the corrections depend 
both on m t , and on while the corrections are proportional to m 2 H . In contrast 
to all previous corrections, they arise due to the electroweak interaction in two loops, 
not one. 

4. In the leading approximation in t the correction S^V^t, h) produced by the diagrams 
of figures 8d-f is independent of i and takes the form 



where the function A(h/t), calculated in [53] , is given in Table 4 for m^/m t < 4. 
For m t = HhGeV and m H = ZOQGeV 

S^Vi = -0.11 . (80) 

The following expansion holds for m^/m ( > 4: 

49 27 3 1 

A(h/t) = it 2 lnr - -ln 2 r - -r(2 - 12tt 2 + 121nr - 271n 2 r) 

v ' ' 4 2 2 3 V ; 



r 



2 



1613 - 240vr 2 - 1500 In r - 720 In 2 r) , (81 



48 



where r = t/h. 5^Vi(t,h) is the greatest of the two-loop corrections in electroweak 
interaction; however, it is also several times smaller than the main gluon corrections 
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Table 4 



m H /m t 


A(m H /m t ) 


r (2) (m H 1 'm t ) 


m H /m t 


A(m H /m t ) 


T {2) (m H /m t ) 


0.00 


0.739 


5.710 


2.10 


9.655 


1.373 


0.10 


1.821 


4.671 


2.20 


9.815 


1.421 


0.20 


2.704 


3.901 


2.30 


9.964 


1.475 


0.30 


3.462 


3.304 


2.40 


10.104 


1.533 


0.40 


4.127 


2.834 


2.50 


10.235 


1.595 


0.50 


4.720 


2.461 


2.60 


10.358 


1.661 


0.60 


5.254 


2.163 


2.70 


10.473 


1.730 


0.70 


5.737 


1.924 


2.80 


10.581 


1.801 


0.80 


6.179 


1.735 


2.90 


10.683 


1.875 


0.90 


6.583 


1.586 


3.00 


10.777 


1.951 


1.00 


6.956 


1.470 


3.10 


10.866 


2.029 


1.10 


7.299 


1.382 


3.20 


10.949 


2.109 


1.20 


7.617 


1.317 


3.30 


11.026 


2.190 


1.30 


7.912 


1.272 


3.40 


11.098 


2.272 


1.40 


8.186 


1.245 


3.50 


11.165 


2.356 


1.50 


8.441 


1.232 


3.60 


11.228 


2.441 


1.60 


8.679 


1.232 


3.70 


11.286 


2.526 


1.70 


8.902 


1.243 


3.80 


11.340 


2.613 


1.80 


9.109 


1.264 


3.90 


11.390 


2.700 


1.90 


9.303 


1.293 


4.00 


11.436 


2.788 


2.00 


9.485 


1.330 









5. Corrections 5 5 Vi due to two-loop diagrams of the type of figure 8g-8j. They are negli- 
gible, but for the sake of completeness of the presentation, we list them in Appendix 
P. 

Accidental (?) compensation and the mass of the t-quark. 

Now that we have expressions for all terms in formula (|58"D, it will be convenient to 
analyze their roles and the general behavior of the functions Vi(t, h). As functions of m t at 
three fixed values of run, they are shown in figures 9a, 10a, 11a. On all these figures, we see 
a cusp at m t = rriz/2. This is a typical threshold singularity that arises when the channel 
Z — > ti is opened. It is of no practical significance since experiments give m t ^> m^/2. 
What really impresses is that the functions V{ are nearly zero in the interval m t ~ 100 — 200 
GeV. This happens because of the compensation of the leading term t and the rest of the 
terms which produce a negative aggregate contribution. This is especially well pronounced 
in the function Vr at m t ~ 180 GeV. Here the main negative contribution comes from the 
light fermions (the constant Cr). 

If we neglect the small correction S^V^t, h) which depends both on t and on h, then each 
function Vi(t,h) is a sum of functions one of which is t-dependent but independent of h, 
while the second is /i-dependent but independent of t (plus, of course, a constant which is 
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50 75 100 125 150 175 200 225 250 

m lop (GcV) 




100 200 300 400 500 600 700 800 900 1000 

m„ (GcV) 



Figure 9: V m as a function of m t for three values of m#: 60, 300 and 1000 GeV, according 
to equation d58|) . The dashed-curve parabola t = {m t /mz) 2 corresponds to the Veltman 



approximation. Solid horizontal line traces the experimental value of V^ p while the long- 
dash horizontal lines give its upper and lower limits of the la level (a). V m as a function of 
ran for three values of m t : 150, 175 and 200 GeV, according to equation flSlp. (Horizontal 
lines the same as in figure 9(a). ) 
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Figure 10: Va as a function of m t . (The remaining clarifications are similar to those given 
to figure 9) (a). Va as a function m# for three values of m t : 150, 175 and 200 GeV. (See 
caption to figure 9) (b) 
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Figure 11: Vr as a function of mt. (The remaining clarifications are similar to those given 
to figure 9) (a). Vr as a function of m# for three values of m t : 150, 175 and 200 GeV. (See 
caption to figure 9) (b). 
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independent of both t and h). Therefore the curves for m# = 60 and 1000 GeV in figures 
9a, 10a, 11a are (mostly) produced by the parallel transfer of the curve for m H = 300 GeV. 

We see in figures 9a, 10a and 11a that if the t-quark were light, radiative corrections 
would be negative, and if it were very heavy, they would be much larger. This looks like a 
conspiracy of the observable mass of the t-quark and all other parameters of the electroweak 
theory, as a result of which the electroweak correction V R becomes anomalously small. Note 
that the corrections do not vanish simultaneously because, if we fix the value of m#, then 
V m (t), V^(t) and Vn(t) cross the horizontal lines V* = at different values of m t . What 
happens is an approximate vanishing of the correction, which as if corresponds to some 
broken symmetry. The nature of this symmetry is not clear at all, and even its existence is 
very problematic. 

One should specially note the dashed parabola in figures 9a, 10a and 11a corresponding 
to the Veltman term t. We see that in the interval < m t < 250 GeV it lies much higher 
than Va and Vr and approaches V m only in the right-hand side of figure 9a. Therefore, the 
so-called non-leading 'small' corrections that are typically replaced with ellipses in standard 
texts, are found to be comparable with the leading term t. 

A glance at figures 9a, 10a, 11a readily explains how the experimental analysis of elec- 
troweak corrections allows, despite their smallness, a prediction, within the framework of 
the minimal standard model, of the t-quark mass. Even when the experimental accuracy of 
LEP1 and SLC experiments was not sufficient for detecting electroweak corrections, it was 
sufficient for establishing the t-quark mass using the points at which the curves Vi{m t ) inter- 
sect the horizontal lines corresponding to the experimental values of Vj and the parallel to 
them thin lines that show the band of one standard deviation. The accuracy in determining 
m t is imposed by the band width and the slope of Vi(m t ) lines. 

The dependence Vi(m#) for three fixed values of m t = 150, 175 and 200 GeV (figures 
9b, 10b, lib) can be presented in a similar manner. As follows from the explicit form of the 
terms Hi(m H ), the dependence Vi(m H ) is considerably less steep (logarithmic). This is the 
reason why the prediction of the higgs mass extracted from electroweak corrections has such 
a high uncertainty. We will see later (figures 13, 14, 15) that the accuracy of prediction of 
m H will greatly depend on what the t-quark's mass is going to be. If m t = 150 ±5 GeV, then 
m H < 200 GeV at the 3a level. If m t = 200 ± 5 GeV, then m H > 120 GeV at the 3a level. 
If, however, m t = 175 ± 5 GeV, we are hugely unlucky: there is practically no constraint on 
m H . 

Before starting a discussion of hadronic decays of the Z boson, let us 'go back to the 
roots' and recall how the equations for Vi(m t ,m H ) were derived. 

How to calculate Vj? 'Five steps'. 

An attentive reader should have already come up with the question: what makes the am- 
plitudes of the lepton decays of the Z boson in the one-loop approximation depend on the 
self-energy of the W boson? Indeed, the loops describing the self-energy of the W boson ap- 
pear in the decay diagrams of the Z boson only beginning with the two-loop approximation. 
The answer to this question is this. We have already emphasized that we find expressions for 
radiative corrections to Z-boson decays in terms of a, mz and G^. However, the expression 
for G M includes the self-energy of the W boson even in the one-loop approximation. The 
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point is thus in our expressing some observables (in this particular case, mw /mz, 9a, 9v / 9a) 
in terms of other, more accurately measured observables (a, mz, G>)- 

Let us trace how this is achieved, step by step. There are altogether 'five steps to 
happiness', based on the one-loop approximation. 

Step I. We begin with the electroweak Lagrangian after it had undergone the sponta- 
neous violation of the SU(2) x [/(l)-symmetry by the higgs vacuum condensate (vacuum 
expectation value - VEV) r\ and the W and Z bosons became massive. Let us consider the 
bare coupling constants (the bare charges e of the photon, g of the VT-boson and f of the 
Z-boson) and the bare masses of the vector bosons: 

mzo = ^foV , (82) 

m W o = ^9oV , (83) 

and also bare masses: m t0 of the t-quark and of the higgs. 

Step II. We express a, G^, m z in terms of / , go, e , rj, m t0 , m m and 1/e (see Appendix 
E). Here 1/e appears because we use the dimensional regularization, calculating Feynman 
integrals in the space of D dimensions (see Appendix A). These integrals diverge at D = 4 
and are finite in the vicinity of D = 4. By definition, 

2e = 4 - D -> . (84) 

Note that in the one-loop approximation m t0 = m t , rrino = %, since we neglect the 
electroweak corrections to the masses of the t-quark and the higgs. For the higgs this 
approximation is quite legitimate, since the accuracy of extracting its mass from radiative 
corrections is very poor. As for the t-quark, this statement is also true at the current 
accuracy of the experimental measurement of electroweak corrections; however, this would 
become an unacceptably crude approximation if this accuracy could be improved by an 
order of magnitude at LEP and SLC The situation here is analogous to that for G M and the 
self-energy of the W^-boson. 

Step II is almost physics: we calculate Feynman diagrams (we say 'almost' to empha- 
size that observables are expressed in terms of nonobservable, 'bare', and generally infinite 
quantities). 

Step III. Let us invert the expressions derived at step II and write fo, go, V in terms of 
a, G^, m z , m t , ran and 1/e. This step is pure algebra. 

Step IV. Let us express V m , Va, Vr (or the electroweak one- loop correction to any other 
electroweak observable, all of them being treated on an equal basis) in terms of fo, go, r], m t , 
mil and 1/e. (Like step II, this step is again almost physics.) 

Step V. Let us express V m , Va, Vr (or any other electroweak correction) in terms of a, 
G M , m z , m t , m H using the results of steps III and IV. Formally this is pure algebra, but in 
fact pure physics, since now we have expressed certain physical observables in terms of other 
observables. If no errors were made on the way, the terms 1/e cancel out. As a result, we 
arrive at formula (|58D which gives Vi as elementary functions of t, h and s. 

The five steps outlined above are very simple and visually clear. We obtain the main 
relations without using the 'heavy artillery' of quantum field theory with its counterterms in 
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the Lagrangian and the renormalization procedure. This simplicity and visual clarity became 
possible owing to the one-loop electroweak approximation (even though this approach to 
renormalizations is possible in multiloop calculations, it becomes more cumbersome than 
standard procedures). As for the QCD-corrections to quark electroweak loops and the two- 
loop higgs contribution hidden in the terms SVi in equation (|58"D , we take the relevant formulas 
from the calculations of other authors. 



Hadronic decays of the Z boson. The leading quarks 
and hadrons. 



As discussed above (see formulas (B5|)-(ED) and the subsequent section on 'Asymmetry'), an 



analysis of hardronic decays reduces to the calculation of decays to pairs of quarks: Z — > qq. 
The key role is played by the concept of leading hadrons that carry away the predominant 
part of the energy. For example, the Z — > cc decay mostly produces two hadron jets flying 
in opposite directions, in one of which the leading hadron is the one containing the c-quark, 
for example, D~ = cd, and in the other the hadron with the c-quark, for example, D° = cu 
or A+ = udc. Likewise, Z — > bb decays are identified by the presence of high-energy B or B 
mesons. If we select only particles with energy close to mx/2, the identification of the initial 
quark channels is unambiguous. The total number of such cases will, however, be small. 
If we take into account as a signal less energetic .B-mesons, we face the problem of their 
origin. Indeed, a pair bb can be created not only directly by a Z boson but also by a virtual 
gluon in, say, a Z —>■ cc decay (figure 12a) or Z — > uu, or ss. This example shows the sort 
of difficulty encountered by experimentalists trying to identify a specific quark-ant iquark 
channel. Furthermore, owing to secondary pairs, the total hadron width is not strictly equal 
to the sum of partial quark widths. 

We remind the reader that for the partial width of the Z — > qq decay we had (j35"D 



T q = T(Z - qq) = 12T [g 2 Aq R Aq + g 2 Vg R 2 Vq ] , (85) 
where the standard width T is, according to (|31|). 

To = = 82.944(6)MeV , (86) 

and the radiators RA q and Ry q are given in Appendix F. As for the electroweak corrections, 
they are included in the coefficients g Aq and gy q . The sum of the Born and one-loop terms 
has the form 

g Aq = T 3q [l + -^V Aq (t,h)] , (87) 

310 I 

R q ee g Vq /g Aq = 1 - 4|Q 9 |s 2 + - 2 * aV Rq (t, h) . (88) 



Air(c 



Decays to pairs of light quarks. 
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Figure 12: The Z — > cc decay producing a secondary pair 66 created by a virtual gluon g 
(a). The Z — > gg decay with a virtual gluon g which connects a side of the quark triangle 
with the external quark line. Diagrams of this type have not been calculated yet (b). The 
verted electroweak diagram involving t-quarks and contributing to the Z — > bb decay (c). 
One of the diagrams describing gluon corrections to the diagram of figure 12(c). 
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Here, as in the case of hadronless observables, the quantities V that characterize correc- 
tions are normalized in the standard way: V — > t as t >> 1. Naturally, those terms in V that 
are due to the self-energies of vector bosons are identical for both leptons and quarks. The 
deviation of the differences V Aq — V A and Vn q — Vr from zero are caused by the differences 
in radiative corrections to vertices Z — > qq and Z — > 11. For four light quarks we have 

1287rs 3 r 3 

V Au (t, h) = V Ac (t, h) = V A (t, h) + — (Fai + F Au ) , (89) 

da 

1287TS 3 C 3 

V Ad (t, h) = V As (t, h) = V A (t, h) + — — (Fai - F Ad ) , (90) 

6a 

^ u (t, fc) = Vkc(*, h) = V R (t, h) + { - 1 x (91) 

da 

x [F vl - (1 - 4s 2 )F Al + |(-(1 - ^s 2 )F Au + i^Vu)] , 

IW* , fc) = Vk.(t, fc) = V R (t, h) + 167lSC{ f_ - ^ x (92) 
x [F vl - (1 - As 2 )F Al + 3((1 - ^s 2 )F Ad - F Vd )\ , 



where (see Appendix M): 



F M = ^(d.0099 + 16.45s 2 ) , (93) 
An 

a 

F vl = — (3.1878 + 14.9<5s 2 ) , (94) 

47T 

F Au = -^(2.6802 + 14.75s 2 ) , (95) 
Air 

F Vu = --^(2.7329 + 14.25s 2 ) , (96) 
An 

a 

F Ad = —(2.2221 + 13.55s 2 ) , (97) 
An 

F Vd = -^(2.2287 + 13.55s 2 ) . (98) 
An 

The values of F are given here for s 2 = 0.23110 — 5s 2 . The accuracy to five decimal places 
is purely arithmetic. The physical uncertainties introduced by neglecting higher-order loops 
manifest themselves already in the third decimal place. It is necessary to point out that 
corrections of the type of that shown in figure 12b have not yet been calculated. 

Decays to bb pair. 

In the Z — > bb decay it is necessary to take into account additional t-dependent vertex 
corrections: 

8s 2 c 2 

V Ab (t, h) = V Ad (t, h) - _ 2g2) (0(t) + 50(t)), (99) 
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V m (t,h) = V Rd (t, h) 



As\c 2 - s 2 ) 
3(3 - 2s 2 ) 



(cj>(t) + 5cj>(t)). 



(100) 



Here the term calculated in |57]] corresponds to figure 12c and the term <50(t) calculated 
in |58| , |53| corresponds to the leading gluon and higgs corrections to the term (f){t) (see 
figure 12d). Expressions for 0(t) and 5<p(t) are given in Appendix N. For m t = 175 GeV, 
d s (m z ) = 0.125, m H = 300 GeV, 

(j)(t) = 29.9 , (101) 



and correction terms in equations 
respectively. 



S(f>(t) = -3.0 , (102) 
|) and ( |100|) are very large: they equal —5.0 and —1.8, 



6 Comparison of theoretical results and experimental 
LEP1 and SLC data. 

LEPTOP code. 

A number of computer programs (codes) were written for comparing high-precision data 
of LEP1 and SLC. The best known of these programs in Europe is ZFITTER |59|] , which takes 
into account not only electroweak radiative corrections but also all purely electromagnetic 
ones, including, among others, the emission of photons by colliding electrons and positrons. 
Some of the first publications in which the t quark mass was predicted on the basis of 
precision measurements pOl, were based on the code ZFITTER. Other European codes, 



BHM [61], WOH, TOPAZO [62], somewhat differ from ZFITTER. The best known in the 



USA are the results generated by the code used by Langacker ] |35|] . 

The original idea of the authors of this review in 1991-1993 was to derive simple analyt- 
ical formulas for electroweak radiative corrections, which would make it possible to predict 
the i-quark mass using no computer codes, just by analyzing experimental data on a sheet 
of paper. Alas, the diversity of hadron decays of Z bosons, depending on the constants of 
strong gluon interaction a s , was such that it was necessary to convert analytical formulas 



into a computer program which we jokingly dubbed LEPTOP [63]. The LEPTOP calcu- 
lates the electroweak observables in the framework of the Minimal Standard Model and 
fits experimental data so as to determine the quantities mt, Mh and a s (Mz)- The logical 
structure of LEPTOP is clear from the preceding sections of this review and is shown in 
the Flowchart. The code of LEPTOP can be downloaded from the Internet home page: 
http: / /cppm.in2p3.fr./leptop/intro_leptop.htmJ 



A comparison of the codes ZFITTER, BHM, WOH, TOPAZO and LEPTOP carried 
out in 1994-95 [|18|] has demonstrated that their predictions for all electroweak observables 



coincide with accuracy that is much better than the accuracy of the experiment. The results 
of processing the experimental data using LEPTOP are shown below. 
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Flowchart. 



Choose three observables measured with the highest accuracy: 

G^Mz ,a{M z ) = a 



Determine angle 9 ( s = sin c = cos 9) 
in terms of G M , M z , a: G M = / y/2)a/ s 2 c 2 M 2 z 



Introduce bare coupling constants in the framework of MSM 

(a , a Z o , (two), bare masses (M Z o, M wo , M HQ , m i0 , m q0 ) 
and the vacuum expectation value (VEV) of the higgs field r\. 



Express ao , azo , Mzo in terms of , Mz , a 
in one-loop approximation, using dimensional 
regularization {1/e 



Express one-loop electroweak corrections to all 
electroweak observables in terms of ao , azo , Mzo , ™>t ? Mff? 
and hence, in terms of G^ , , a ,m t , Mjj. Check 
cancellation of the terms , /i). 



Introduce gluon corrections to quark loops 
and QED (and QCD) final state interactions 
(in hadron decays), given in terms of 
a ,a s (M z ) ,m b (M z ) ,m t . 



Compare the predictions of the Born approximation 
and the approximation Born + one loop with the experimental 

data on Z-decays. 



Perform global fit for the three parameters 
mt , &s(Mz) , Mh or for the first two, with the third fixed. 



Derive theoretical predictions of the central values 

for all electroweak observables and of 
the corresponding uncertainties ('errors'). 
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General fit. 



Table 5 shows experimental values of the electroweak observables, obtained by averaging 
the results of four LEP detectors (part a), and also SLC data (part b) and the data on W 
boson mass (part c). (The data on the W boson mass from the pp-colliders are also shown, 
for the reader's convenience, in the form of s^y, while the data on s 2 ^ from ^iV-experiments 
are also shown in the form of mw- These two numbers are given in italics, emphasizing that 
they are not independent experimental data.) We take experimental data from the paper 

Ill- 
Table 5 sums up the experimental data used for determining (fitting) the parameters of 

the standard model m t and a s {m z ) (see Table 7). The central values in the third column of 

Table 5 were calculated for = 300 GeV. Shown in brackets is the uncertainty of the last 

significant decimal places due to the uncertainty of the fitted values m t and a s . Above and 

below we give the shifts in the last significant decimal places corresponding to itlh = 1000 

GeV and rau = 60 GeV, respectively The last column shows the value of the 'pull'. By 

definition, the pull is the difference between the experimental and the theoretical values 

divided by experimental uncertainty. The pull values show that the maximum discrepancy 

between the experimental data and the Minimal Standard Model is found for (3.8a). 

Deviations at the level of 2.5a are also found for R c and sin 2 O^fl = sf in A^r (SLC). For 

most observables the discrepancy is less than la. At the same time, Table 6 shows that the 

value sf = 0.23186(34) extracted from all asymmetries at LEP agrees quite well with the 

fitted MSM value sf = 0.2321(4) from the LEP data, and for all sets of data. 
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Table 5 



Observable 


Experimental 
data 


Fit 

standard 
model 


Pull 


a) LEP 








shape of Z-peak 

and lepton asymmetries: 

m z [GeV] 

T z [GeV] 

a h [nb] 

Ri 

A 1 

FB 

r-polarization: 

A e 

Results for b and c 

quarks: 

Rb 

R c 

A b FB 

A c 

FB 

Charge asymmetry for pairs 

£ 1 * 1 j_ 1 — 

or light quarks qq: 
sf{< Qfb >) 


91.1884(22) 
2.4963(32) 
41.488(78) 
20.788(32) 
0.0172(12) 

0.1418(75) 
0.1390(89) 

0.2219(17) 
0.1543(74) 
0.0999(31) 
0.0725(58) 

0.2325(13) 


2.4976(26)^ 6 
41.450(20);? 
20.770(24)+^ 
0.0158(6);^ 

0.1450(26);L3 
0.1450(26);? 3 

0.2155(3);? 

0.1724(1)11 
0.1017(19);? 
0.0726(14);^ 

0.2318(3) + i 


-0.4 
0.5 
0.6 
1.2 

-0.4 
-0.7 

3.8 
-2.5 
-0.6 
0.0 

0.6 


b) SLC 

A LR 

sf(A LR ) 

Rb 
A b 
A c 


0.1551(40) 
0.2305(5) 
0.2171(54) 
0.8410(530) 
0.6060(900) 


0.1450(26);I 3 
0.2318(3)±\ 
0.2155(3);? 
0.9345(2J + 4 

0.6670(12);^ 


2.5 

-2.5 

0.3 

-1.8 

-0.7 


c) pp and uN 








m w [GeV] (pp) 
\-m 2 w jm\ (uN) 


80.26(16) 
0.2253(31) 
0.2257(47) 
80.24(24) 


80.35(5)^ 
0.2237(9);^ 


-0.5 
0.4 
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Table 6 



Observable 


2 


Average over 
groups of observations 


Cumulative x /d.o.j. 
average 


A 1 

^-FB 

A T 
A e 


0.23096(68) 
0.23218(95) 
0.2325(11) 


0.23160(49) 


0.23160(49) 1.9/2 


^FB 

A c 
FB 


0.23209(55) 
0.2318(13) 


0.23205(51) 


0.23182(35) 2.4/4 


< Qfb > 


0.2325(13) 


0.2325(13) 


0.23186(34) 2.6/5 


A LR (SLD) 


0.23049(50) 


0.23049(50) 


0.23143(28) 7.8/6 



Table 6 gives experimental values of sf. The third column was obtained by averaging 
of the second column, and the fourth by cumulative averaging of the third; it also lists the 
values of x 2 over the number of degrees of freedom. 



Table 7 



Physical 
quantities 


LEP 


LEP + SLC 


LEP + M w 


LEP + SLC 

+ M w 


m t (GeV) 


171(9)i« 21 


182(7)31 


170(8)3! 


180(7)t£ 


a s 


0.125(4)3 


0.123(4)3 


0.125(4)3 


0.124(4)1^ 


X 2 /d.o.f. 


18/9 


29/13 


18/11 


30/15 


s 2 


0.2321(4)3 


0.2317(3)3 


0.2321(4)3 


0.2318(3)3 


gvi/gAi 


0.07116(16)+| 


0.0732(12)+! 


0.0716(16)+! 


0.0728(12)+! 


s 2 


0.2247(9)3 


0.2234(9)3 


0.2237(9)3 


0.2237(9)3 


m w /m z 


0.8805(5)3 


0.8813(5)3 


0.8804(5)3 


0.8811(5)3 


m w (GeV) 


80.29(5)3 


80.36(5)3 


80.28(5)3 


80.35(5)3 



Table 7 gives the fitted values m t and a s = a s (m 2 z ), and also the values of x 2 times the 
number of degrees of freedom for various sets of data, where M w stands for both the data 
of M\y measurements in pp collisions and the values of s 2 ^, extracted from experiments with 
vN. The lower part of the table gives the values of sf = sin 2 O^ff = |(1 — gvi/dAi) and 
sfy = 1 — m^/ml, calculated in one-loop electroweak approximation in the framework of 
the Minimal Standard Model using fitted values of m t and a s . Errors given in parentheses 
are due to uncertainties in m t , a s and a, and were calculated by summation of squares, 
ignoring correlations. Note that the errors in the values of s 2 ^ calculated using the fits are 
substantially lower than in the experimental values of this quantity (see Table 5); at the 
same time, the errors for sf are practically identical for the experimental (Table 6) and the 
theoretical (Table 7) values. Note that the first and second rows of the lower part of Table 
7 carry identical information; the same is true for the third, fourth and fifth rows. 
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7 Conclusions. 



Achievements. 

What are the main results obtained with four detectors of the LEP1 collider and one 
detector of the SLC? 

Judging by the criteria of accelerator and experimental techniques, the highest imaginable 
level has been achieved in the team creativity. The impossible became the reality owing to 
a never before dreamt-of sophistication of equipment of gigantic high-precision detectors. 
Twenty million decays of Z bosons were measured with better accuracy than that admired 
in gemstone cutting. 

From the physics standpoint, the main result is the experimental proof that there exist 
only three generations of quarks and leptons with light neutrinos. The number of light 
neutrinos, found from the sum width of the invisible Z boson decays is 

N v = 2.990(16) (103) 

The lower limit on the masses of the heavy neutrinos in the additional generations, provided 
they exist, is close to m z j1 and equals 44 GeV. 

No new particles were found in Z boson decays. In particular, the higgs was not found. 
From LEP1 data, the lower limit on the higgs mass is 

m H > 60 GeV . (104) 

The high-precision measurement of the Z boson mass, its total and partial decay widths 
and also of the P- and C-violating asymmetries made it possible to experimentally determine 
the electroweak radiative corrections. A comparison of these experimental values with the 
results of theoretical calculations based on the Minimal Standard Model led to prediction of 
the t-quark mass m t and the constant of strong interaction for gluons a s \ 

m t = 180(7)t^GeV , (105) 



a s = 0.124(4)^ . (106) 

Shown in parentheses here is the uncertainty (one standard deviation) due to the uncertainty 
of the experimental data. The central value corresponds to the assumption that tuh = 
300 GeV, and the upper and lower 'shifts' correspond to m H = 1000 GeV and 60 GeV, 
respectively. Radiative corrections depend only weakly on m#, so one cannot extract the 
value of m,H from them. The t-quark mass predicted on the basis of the radiative corrections 
within the currently known experimental uncertainties is in excellent agreement with the 
results of direct measurements of m t in CDF and DO detectors at the Tevatron []: 

m t = 176(13) GeV [13] CDF . (107) 
1 At the spring 1996 conferences more accurate data have been presented: 

m t = 175.6 ± 5.7 ± 7.4 GeV CDF@ 
mt = 170 ± 15 ± 10 GeV D0|| 
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m t = 199(30) GeV [14] DO . (108) 

After the t-quark mass uncertainty is further reduced, radiative corrections can be used 
for determining the region in which the higgs mass can lie. Figures 13-15 show that the 
results will greatly depend on luck. If m t = 150(5) GeV, figure 13 demonstrates that 
m H < 150 GeV at the 3a level. If m t = 200(5) GeV, then m H > 120 GeV at the 3a level 
(figure 15). If, however, m t = 175(5), then the higgs can have any mass within 3a. 

Problems. 

A cursory glance at Table 5 is sufficient for identifying the main problem of the Z- 
boson physics: the discrepancy between the measured width of decay into a pair bb and its 
theoretical prediction. 

The supersymmetrization of the standard model may help solving this problem |66[ (see 
figure 16a, this diagram with superpartners increases Rb ). Theories with extra neutral 
vector bosonZ also may be relevant. Let us turn now to a s (mz)- A number of papers 
67] , |68| pointed out that the value 0.124(4)12, shown in Table 7, is in contradiction with 



the measurements of a s (q 2 ) for q 2 < (lOGeV) 2 in deep inelastic scattering [|69[] , in hadron 



decays of the T-meson |7(| and especially in the spectrum of upsilonium levels [[7IJ. If the 
low-energy values a s (q 2 < (lOGeV) 2 ) are extrapolated in the framework of the standard 
QCD to q 2 = m|, we find a s (m 2 z ) = 0.110 — 0.118. As for the uncertainty in this range, 



the authors of p8]|-||71| do not come to a common opinion. The most cautious evaluate it 



as ±0.005 US9[|-[f70"l|. The bravest one insists on ±0.001 f68[, flnH . In the last case, there is 



an obvious contradiction with the value derived by analyzing the Z boson decay data. This 
contradiction served as a basis for hypotheses [^7] that MSM predictions for hadron widths 



are modified by the contribution of some new unknown particles to electroweak radiative 
corrections, for example, squarks and gluino, that is, the supersymmetric partners of quarks 
and gluons. For loops with these particles to result in sufficiently strong deviations from 
MSM, it is necessary that squarks and gluino were sufficiently light, with masses of the 
order of 100 GeV. A deviation of the observable value of Rb by more than 3a from the 
values predicted in MSM on the basis of the global fit (see Table 5) is also considered as an 
independent argument in favor of the view that the new physics 'lurks round the corner'. 
Two other deviations from MSM in Table 5 are less serious: R c and A^r are off by 2.5 
standard deviations. Note that in the latter case we witness a discrepancy not only with 
MSM but also with the measurements of the lepton asymmetries at LEP1, since, according 
to equation (|4"7j), A LR = A\. 

Various parametrizations of the manifestations of the new physics can be found in the 
literature. The better known ones are the parameters S,T, U |72[ and Ei, e 2 , £3, £b |73|| . 



Prospects. 

The LEP collider completed its work in the LEP1 mode in autumn 1995 and began 
working in the LEP2 mode. It is expected that the total energy of the electron and positron 
collision will be raised to 192 GeV. What are the main goals of LEP2 [74] ? 

A careful measurement of the cross section of the reaction e + e~ —>■ W + W~ needed to 
measure the P^-boson mass with accuracy of the order of 50 MeV and to test whether the 
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Figure 13: Isolines of x 2 m the m t ,m h plane, obtained by fitting the electroweak corrections 
under the assumption that direct measurements of the t quark mass will give m t = 150 ± 5 
GeV . 
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Figure 14: Same as in figure 13, for m t = 175 ± 5 GeV. 
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Figure 15: Same as in figure 13, for m t = 200 ± 5 GeV. 
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interaction of W bosons with photons and with the Z boson agrees with the standard model 
(see figure 16b). 

A search for a higgs with a mass up to 92 GeV in the reaction e + e~ — > HZ (figure 16c). 
A search for light superparticles (sleptons, squarks). 
A search for the unanticipated. 

Further prospects for testing the Standard Model and for searching for 'new physics' 
beyond its limits are rooted in the Large Hadron Collider (LHC) (a decision to build it has 
already been made at CERN) and in the so-called 'Next Linear Collider' for electrons and 
positrons, which is so far at the stage of discussion of competing projects. 

We are grateful to A. V. Novikov for help in preparing this review; A. R. is grateful 
to CNRS-IN 2 P3/CPPM for support; M.V., V.N. and L.O. are grateful to RFBR for the 
grant 93-02-14431; M.V. and V.N. are grateful to INTAS for the grants 93-3316 and 94-2352, 
which helped to carry out the work presented here. 
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Appendix A. 
Feynman rules in electroweak theory. 



A consistent derivation of Feynman rules for theories with spontaneous violation of gauge 
symmetry can be found in textbooks (see, for example, Itzykson, Zuber; Ramond; Slavnov, 
Faddeev |3[]). In this appendix we only give a summary of results, accompanied by brief 
comments. 

Gauges and propagators. 

Quantization of gauge fields (in MSM this is W^, and A^) requires fixing gauge. The 
most popular is the so-called gauge which corresponds to adding gauge fixing new terms 
SCqf to the classical Lagrangian: 



(A.l) 



$ " I -j,(rj + H(x) + iG°(x)) ) (A ' 2) 



where G ± , G° and H are the components of the higgs doublet $ in the parametrization 

G+(x 

In what follows we use the particular case of gauge, namely 

t,A = = iz = £ 

With gauge fixed, it is possible to determine the propagators Dj^(p), D^ u (p) and D^ v {p) of 
the fields W^, and A^. 



p 2 -my V + ie { p 2 -mf v ^ + ie 



P^Pu, -P^Pu 1 



(A.3) 



[9^ ~ -75-) - z 



9 2 1 1 ) 2 9 2^1- 

jrr - 77% + «e 77% 77% p z - r?%4 + «e 



Dl{p) 



p A 



m z + ie 



p 2 — m z + ie 



(9, 



9iiv 



ml ' 



;i-o 



PyPv 

p 2 — m z £ + ie 



(A.4) 



m z P 2 ~ m l£ + i> e 
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^M = ~{^-^-i)f^-} , (A.5) 

The case £ = 1 corresponds to the Feynman-'tHooft gauge, £ = to the Landau gauge, 
£ — > oo to the Proca gauge. 

As follows from (A. 3) and (A. 4), the propagators of massive vector fields can be written 
as sums of a propagator in the Proca gauge that describes the propagation of physical degrees 
of freedom of a vector particle, and a scalar propagator with a gauge- dependent pole which 
corresponds to the propagation of non-physical degrees of freedom. As a result, diagrams 
with virtual W ± , Z bosons contain non-physical threshold singularities whose positions 
depend on the gauge parameter £. These non-physical singularities partially cancel out after 
the appropriate diagrams with virtual Goldstone bosons G ± , G° (arising from the higgs 
doublet $ (|A.2|) ) are added. The Goldstone boson propagators have the form: 

D G+ (P) = - ( A - 6 ) 

p z - r7%£ + is 

D G°(P) = ~ % TT-~- (A-7) 

p z — m z z t, + %e 

A complete restoration of unitarity (cancellation of non-physical singularities) and of 
gauge invariance (validity of the Ward identities) are achieved if one takes into account the 
diagrams with the Faddeev-Popov ghosts r] ± , r] z and r] A , which interact only with gauge 
fields and the Goldstone fields and which do not correspond to any physical degrees of 
freedom. 

Ghost propagators take the form 

= p2 _4; + fe ' (A - 8) 

A,«(P) = (A. 10) 

p z + IE 

Ghosts obey the Fermi statistics, so an additional sign (—1) must be ascribed to ghost loops, 
as one does for fermion loops. 

The propagators of other fields are written as 

for the higgs field : Duip) = n (A. 11) 

p z — m H + ie 

for fermion fields : Df(p) = (A. 12) 

p — rrif + ie 

In order to describe numerous three-particle vertices, it is convenient to unify the notations. 
Let us fix the momenta once and for all, as shown in figure 17, and let us denote a vertex by 
a set of fields in the following order: (ACB). The Feynman rules for three-particle vertices 
are then written as follows. 
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<^> (ABC) 



B 



Figure 17: Three-particle vertex. 



Interaction between gauge fields and fermions. 



(/A*/) 



(fZ»f) -i{\9vlp + 9AlvTl!>\ , 



(A.13) 



(UW-D) -i^=V DU ^(l + 15 ) ; 

where Qf is the charge of the fermion /, T[ is the third projection of the isotopic spinor, 
describing the left-handed component of the fermion /, = 2T 3 ^, gy = 2Tl — 4Q/sin 2 #; 
U and D denote any of the quarks with T/ = \ and T( = — |, respectively, and Vdu is an 
element of the Kobayashi-Maskawa matrix. 

Interaction of scalar fields with fermions. 



(fHf) 
(fG°f) 
(UG~D) 
{DG+U) 



2mw 



9 m i rpf 
4 mw ^ ' 



-^J^[( m D ~ mu) + l5 (m D + mu )} , 
~172^( mD ~ miJ ) -ls( m D +mu)} . 



(A.14) 
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Three-boson interactions. 

Three gauge bosons: 

(WfA v W+) ie[r + q) x g^ - (q + p) v g Xii + (p - r)^g vX \ , 

(A.15) 

(W x ^Z l/ W+) ig cos 0[r + q)xg llv - (q + p)„^ + (p - r)^g uX ] 
Two gauge bosons and one scalar boson (G or H): 

(W+G~A V ) iem w g^ v , 

{W+G~Z U ) -igm z sin 2 9 g^ , 

(A.16) 

(W+HW+) 

2 

(Z^HZ U ) ig^g^u ; 

One gauge boson and two scalar bosons (GG or GH): 



(G+A.G+) 


-ie(p + q)^ , 


(G+Z.G+) 


-Wf^(p + q)» 


{G°W+G+) 


\g(p + q)n , 


(HW+G+) 


-\ig{p + q)^ , 


(G+W^G ) 


~h(p + q)n , 


(G + W~H) 


-\ig{p + q)» , 


(HZ^G ) 


^ke(p + q)n ; 



Interaction of higgses with goldstones and among themselves: 

(G+HG+) ig& , 

(G-HG-) -\g^ , 

(A.18) 

(G°HG°) -ig& , 
(HHH) -%g&L . 



53 



Interaction between ghosts and gauge fields: 

{^A^) -ieq^ , 

{r]-W+rfi) -ieq^ , 

(ri + W~ri~<) ieq^ , 

{rfW-rf) -ieq^ , 

(ri + Z^ + ) -ig cos 9 q^ 

(rj-Z^rj-) ig cos 9 q^ , 
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<^> (ABCD) 



Figure 18: Four-particle vertex. 

(V~W+T] Z ) -ig cos9 , 

(r) + W-r) z ) ig cos 9 q^ , 

{rfW+rf) ig cos 9 q^ , 

{rfW'Tf) -ig cos 9 q^ ; 
Interaction of ghosts with a higgs or a goldstone: 

(rj-Hrj-) -\gim w , 

(r]-G°r]-) -\gim w , 

( V + G° V + ) \gim w , 

{rfG + rj + ) —ie^mw , 

(V Z G + T) + ) -- 2 9 C ^rn w , 

(r]-G + r] z ) \gim z . 



(A.20) 



Four-boson interactions. 



To describe four-boson vertices, we introduce the notation (ABCD), see figure 18. In this 



55 



(A.21) 



notation, the interactions of four vector bosons take the form 

(W^W-A^Ap) -ie 2 [2g upgilX - g w g vX - g^gxp] , 

( w x +W p Z uZ P ) -ig 2 cos 2 9[2g up g flX - g w g uX - gpu9x P ] , 

( w x +W p A uZ p ) -ieg cos O^g^g^ ~~ 9 PP 9vX - 9^9x P ] ■ 
Interaction between two vector bosons and HH, GG or HG. 

(W+W-HH) l -g 2 g, u , 

{W;W;G"G") l -g 2 g, u , 

(W+W;G-G + ) % -g 2 g, v , 

(A,A V G-G + ) 2ie 2 g, v . 

Z^Z V HH i^g 2 sec 2 9g pu , 

Z,Z U G'G° t^g 2 sec 2 9g pu , 

Z^Z„G-G + i l -g 2 sec 2 9 cos 2 29g pu , (A.22) 

A^W+G+H i l -egg pv , 



A,W+G + G° -\egg, u , 

A^W-G-G° X -egg, v , 

A P Z U G~G + ieg sec 9 cos 26 ] g^ v , 

Z^W+G+H i^g 2 sec9(^cos29-l)g^ , 

Z^W-G-H i^g 2 sec0(^cos20-l)g^ , 

Z^W+G + G° -^g 2 sec9^cos29-l)g^ , 



Z^W-G-G° l -g 2 sec 0(± cos 29 - l) 9ilv 
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Interactions of GGGG, HHHH or GGHH. 

(G + G-G-G + ) -y 2 w? H /m 2 w , 
-{g 2 m 2 H /m 2 w , 
(G+G-HH) -\g 2 m 2 H /m 2 w , 

(A.23) 

(HHHH) -i\9 2 m 2 H /m 2 w , 

(G G°G°G ) -i\g 2 m 2 H /m 2 w , 
(G°G°HH) -{g 2 m 2 H /m 2 w . 
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Regularization of Feynman integrals. 

Integrals corresponding to diagrams with loops formally diverge and thus need regular- 
ization. Note that there does not exist yet a consistent regularization of electroweak theory 
in all loops. A dimensional regularization can be used in the first several loops; this corre- 
sponds to a transition to a D-dimensional spacetime in which the following finite expression 
is assigned to the diverging integrals: 



d D P ( P 2 y ttt r(f + s )r(a-f -s) 



X 



• (A.24) 



where /x is a parameter with mass dimension, introduced to conserve the dimension of the 
original integral. 

This formula holds in the range of convergence of the integral. In the range of divergence, 
a formal expression ( A.24j) is interpreted as the analytical continuation. Obviously, the 



integral allows a shift in integration variable in the convergence range as well. Therefore, a 
shift p —>■ p + q for arbitrary D can also be done in ( |A.24j ). This factor is decisive in proving 
the gauge invariance of dimensional regularization. 



At D = 4 the integrals in ( |A.24j ) contain a pole term 



2 m 2 

A = - - + hi47r- 7 -ln— , (A.25) 

4 — D /j, z 



where 7 = 0.577... is the Euler constant. Choice of constant terms in (A.25) is a matter of 
convention. 

The algebra of 7-matrices in the D-dimensional space is defined by the relations 

lulu + Ivlp, = 2g IJV x / , (A. 26) 

9w = D , (A.27) 
7„7^ = (2 - D) lv , (A.28) 

where I is the identity matrix. 

As for the dimensionality of spinors, different approaches can be chosen in the continua- 
tion to the D-dimensional space. One possibility is to assume that the 7 matrices are 4x4 
matrices, so that 

SpI = A . (A.29) 

The D-dimensional regularization creates difficulties when one has to define the abso- 
lutely antisymmetric tensor and (or) 75 matrix. For calculations in several first loops, a 
formal definition of 75, 

757^ + 7/.7S = , (A.30) 
75 = I (A.31) 

does not lead to contradictions. 
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Thus, the amplitudes of physical processes, once they are expressed in terms of bare 
charges and bare masses, contain pole terms ~ 1/(D — A). 

If we eliminate bare quantities and express some physical observables in terms of other 
physical observables, then all pole terms cancel out. The general property of renormalizabil- 



ity guarantees this cancellation. (We have verified this cancellation directly in [42|.) This 
renormalization procedure is employed in this review. 

In order to avoid divergences in intermediate expressions, one can agree to subtract 
from each Feynman integral the pole terms ~ 1/(4 — D), since they will cancel out anyway 
in the final expressions. Depending on which constant terms (in addition to pole terms) 
are subtracted from the diagrams, different subtraction schemes arise: the MS scheme 
corresponds to subtracting the universal combination 

2 

7 + In Air . 

A-D 1 
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Appendix B. 
Relation between ol and a(0). 



We begin with the following famous relation of quantum electrodynamics ||75|| : 



Here the fine structure constant a = a(0) is a physical quantity. It can be measured as a 
residue of the Coulomb pole 1/q 2 in the scattering amplitude of charged particles. As for the 
running coupling constant a(q 2 ), it can be measured from the scattering of particles with 
large masses m at low momentum transfer: m 3> J \q 2 \. In the standard model we have the 
Z-boson, and the contribution of the photon cannot be identified unambiguously if q 2 ^ 0. 
Therefore, the definition of the running constant a(q 2 ) becomes dependent on convention 
and on details of calculations. 

At q 2 = m 2 z , the contribution of IV-bosons to a = a(m 2 z ) is not large, so it is convenient 
to make use of the definition accepted in QED: 

where 

6a = -n 7 (m|) + s;(o) , 

n 7 W) = 4-S 7 (m|) . (B.3) 

Tib 7 



One-loop expression for the self-energy of the photon can be rewritten as 76 



S 7 (s) = {a/3ir)J2N!Q 2 f [sA f + {s + 2m 2 f )F{s,m f ,m f ) -3/2] - 
/ 

- (a/47r)[3sA^ + (3s + 4m 2 w ) F (s , m w , m w )} , (B.4) 

where s = q 2 , the subscript / denotes fermions, the sum runs through lepton and quark 
flavors, and N[ is the number of colors. The contribution of fermions to £ 7 (g 2 ) is independent 
of gauge. The last term in ( |B.4| ) refers to the gauge-dependent contribution of VT-bosons; 
the 't Hooft-Feynman gauge was used in equation ( |B.4|) . 
The singular term Aj is: 

Aj = 7 + In 47r — In — |- , (B.5) 

e fi z 

where 2e = 4 — D (D is the variable dimension of spacetime, e — > 0), 7 = — r'(l) = 0.577... 
is the Euler constant and fi is an arbitrary parameter. Both 1/e and \x vanish in relations 
between observables. 
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The function F(s, mi, m 2 ) is denned by the contribution to self-energy of a scalar particle 
at q 2 = s, owing to a loop with two scalar particles (with masses mi and m 2 ) and with the 
coupling constant equal to unity: 



F(s,m 1 ,m 2 ) 



1 + 



mr — m 



2 , m i 
2 log — 

m 2 



dxlog 



x 2 s — x(s + ml — m 2 ,) + m\ — it 



m\m 2 



(B.6) 



The function F is normalized in such a way that it vanishes at q 2 = 0, which corresponds to 
subtracting the self-energy at q 2 = 0: 



F(0,m 1 ,m 2 ) = , 
The following formula holds for mi = m 2 = m: 



(B.7) 



F(s,m,m) = F(t) = 



1 — ^4r — 1 arcsin -7= 



1- v^47ln 1+ ^^ 



4r > 1 , 
4r < 1 , 



where r = m 2 / s. 

To calculate the contributions of light fermions, the t-quark and the W^-boson to 5a, we 
need the asymptotics F(t) for small and large r: 



F(t) ~mr + 2 + ..., |r| < 1 , 
F ( T ) - T~ + 777^ + l r l > 1 , 



F'(s,m,m) = —F(s,m,m) s ~° — 



As a result we obtain 

n 7 (m|) 



6r 60r 2 
S 7 (m|) 



1 1 

+ 



6 30r 



(B.8) 
(B.9) 
(B.10) 



m\ 
a ^2 



7T 

a 

47T 



A t + (1 + 2t)F(t) 



1 



[3A W + (3 + 4c 2 )F(c 2 



where t = m 2 /m 2 z , and 



4:71 



(B.ll) 



(B.12) 



5a 




1 + 2t)F(t) - - 




(B.13) 
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Therefore, Sa is found as a sum of four terms, 



Sa = 5on + 5ah + 5a t + Saw 



(B.14) 



Sai 



Sa t 



a 
3^ 



E 



3 L 



In 



m 



mt 



-5/3 



0.03141 



a 4 



fm z \'- 



-0.00005(1) 



(B.15) 



(B.16) 



where we have used that m t 



n 45 \M t J 

175 ± 10 GeV. Note that Sa t is negligible and has the 



antiscreening sign (the screening of the t-quark loops in QED begins at q ^> mf, while in 
our case q 2 = m| < m 2 ). 
Finally, the W loop gives 



Sa 



w 



a 
2^ 



(3 + 4c 2 ) 1 - V4c 2 - 1 arcsin — - 1/3 



2c 



0.00050 



(B.17) 



The value of Saw depends on gauge JFTjj; here we give the result of calculations in the 
't Hooft-Feynman gauge. Traditionally, the definition of a takes into account the contribu- 
tions of leptons and five light quarks; the terms Sat arid Saw are taken into account in the 
electroweak radiative corrections. In our approach, these terms give the corrections 8{Vi. 
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Appendix C 
Summary of the results for a. 

Among the three input parameters a, and Mz, the first one has the maximum un- 
certainty; this uncertainty leads to the uncertainty ±5 GeV in the value of the t-quark 
mass extracted from the measurements of mw and the decay parameters of the Z-boson. 
According to the formulas of Appendix B, 

a = — , 5a = 5ai + bah , 

1 — oa 

5ai = 0.0314 . 

In (pm, Burkh) 5ah was calculated by substituting experimental data on (Te+e-^hadrons 
into the dispersion integral at ^/s < 40 GeV and the parton model result at y/s > 40 GeV: 

5a [ h Burkh] = 0.0282(9) , a [Burkh] = [128.87(12)]- 1 . (C.l) 



In (|26||, Vys) it was pointed out that the simplest model (the lightest vector meson + 



the QCD-improved parton continuum in each flavor channel) produces a surprisingly close 
result: 

Sa[ Yys] = 0.0280(4) , « [Vys] = [128.90(e)]- 1 . (C.2) 

The model with infinite number of poles ( ||26|| , Gesh) yields the determination of 5at 
with very high accuracy: 

Sa [ h Gesh] = 0.0275(2) , a [Gcsh] = [128.96(3)]- 1 . (C.3) 

A recent analysis of experimental data ( [p6j , Sw, f26fl , Ma) yielded considerably lower 
values of 5ah- 

Sa [ h Sw] = 0.0265(8) , a [Sw] = [129. 10(12)]-! , (C.4) 
5a[ m = 0.0273(4) , a [Ma] = [128.99(e)]- 1 . (C.5) 



In this review we make use of the results of a recent analysis ( ||26|| , Eid), 

Saf id] = 0.0280(7) , a [Eid] = [128.896(90)] . (C.6) 
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Appendix D. 
How OLw{q 2 ) and OLz{q 2 ) 'crawl'. 

The effect of 'running' of electromagnetic coupling constants a(q 2 ) (logarithmic depen- 
dence of the effective charge on momentum transfer q 2 ) is known for more than four decades) 
75|j . In contrast to ct(q 2 ), the effective constants of W and Z bosons aw(q 2 ) and az(q 2 ) in 
the region < q 2 < m| 'crawl' rather than run [7S 



If we define the effective gauge coupling constant g 2 {q 2 ) in terms of the bare charge g$ 
and the bare mass m , and sum up the geometric series with the self-energy S(g 2 ) inserted 
in the gauge boson propagator, this gives the expression 



i/0 q 2 —m 2 

Here m is the physical mass, and S(g 2 ) contains the contribution of fermions only, since 
loops with W, Z and H bosons do not contain large logarithms in the region \q 2 \ < m 2 z . 

The bare coupling constant in the difference l/g 2 (q 2 ) — l/g 2 (0) is eliminated, which gives 
a finite expression. The result is 

l/az(q 2 ) ~ l/az(0) = b z F(x) , where x = q 2 /m 2 z , (D.2) 

V%(? 2 ) - 1/W(0) = b w F(y) , where y = q 2 /m 2 z , (D.3) 

F(x) = — — ln\x\ CD A) 

1 — x 



If x 3> 1, equations ( p.2|) and (p.3| ) define the logarithmic running of charges owing 



to leptons and quarks, and bz and bw represent the contribution of fermions to the first 
coefficient of the Gell- Mann-Low function: 

bz = ^{iV u 3[l + (1 - jjs 2 ) 2 ] + Ar d 3[l + (-l + ^ 2 ) 2 ]+iV,[2 + (l + (l-45 2 ) 2 )]} , 



bw = -^{6^ + 2^] , (D.5) 

l07T 



where N u ^ q> i are the numbers of quarks and leptons with masses that are considerably lower 
than y/q 2 . 

For q 2 < m%, the numerical values of the coefficients bz,w a re [^8 



bz ^ 0.195 



b w ^ 0.239 



The massive propagator q2 } m 2 in (p.l|) greatly suppresses the running of aw(q 2 ) and az(q 2 )- 
Thus, according to ( p.2| ) and ( |D.3| ), the constant az(q 2 ) grows by 0.85% from q 2 = to 



q 2 = m z , 



l/a z {m 2 z ) = 22.905 
l/a z {m%) - l/az(0) = -0.195 , (D.6) 
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and the constant aw{<f) grows by 0.95%, 

l/a w (m 2 z ) = 28.74 

l/a w (m 2 z ) - l/a w (0) = -0.272 , (D.7) 
while the electromagnetic constant a(q 2 ) increases by 6.34%: 

l/a(m 2 z ) - l/a w (0) = 128.90 - 137.04 = -8.14 (D.8) 

With the accuracy indicated above, we can thus assume 

a z (m 2 z ) ~ a z (0) 
a w {m 2 z ) ~ a w (0). (D.9) 

At the same time, a(m 2 z ) differs greatly from a(0); therefore the latter has no connection 
to the electroweak physics but only to the purely electromagnetic physics. 
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Appendix E 

Relation between a, rrtz and the bare quantities. 

The bare quantities are marked by the subscript '0'. In the electroweak theory, three 
bare charges e , fo and g that describe the interactions of 7, Z and W are related by a 
single constraint: 

(e /go) 2 + (9o/ fo) 2 = 1 • (E.l) 

The bare masses of the vector bosons are defined by the bare vacuum expectation value 
of the higgs field r] : 

m Z o = -^foVo , m wo = -g r] . (E.2) 
The fine structure constant a = e 2 /4ir is related to he bare charge eo by the formula 

«-«(g 2 = 0) = |ifl-S' 7 (0)-2^^4^) , (E-3) 
4tt \ ' c m z J 

where S'(0) = lim g 2^ S(g 2 )/g 2 . In the Feynman gauge £ 7 z(0) ~ — (a/2%)(rn 2 v /cs)(l/e), 
where the dimension of spacetime is D = 4 — 2s. In the unitary gauge S 7 ^(0) = 0. 

The simplest way to verify the presence the term 2(s/c)£ T z(0)/m| is by considering 
the interaction of a photon with the right-handed electron e#. Note that in this case there 
are no weak vertex corrections due to the P^-boson exchange. (Note also that the left- 
handed neutrino remains neutral even when loop corrections are taken into account, since 
the diagram with the 7 — Z — vip^ interaction is compensated for by the vertex diagram 
with the W exchange). 

Our first basic equation is the renormalization group improved relation between a = 
a(q 2 = m|) and «o: 

1 - n 7 (m|) - 2-n 7Z (o) 



a = ao 



c 



(E.4) 



where II 7 (g 2 ) = S 7 (g 2 )/m|, II 7Z (g 2 ) = S 7Z (g 2 )/m 



The second basic equation is: 



m 7 = m 



| Q [1 - n z (m|)] = m 2 W0 /cl[l - Uz(m 2 z )] . (E.5) 



A similar equation holds for m 



2 . 

w 



m 2 w = m 2 WQ [l - Ii w (m 2 w )} , (E.6) 

where U^q 2 ) = E^/mf, i = W,Z. 
Finally, the third basic equation is 

G,= A °\ [1 + IM0) + D] , (E.7) 

where 11^(0) = ^w^/m^ comes from the propagator of W, while D is the contribution of 
the box and the vertex diagrams (minus the electromagnetic corrections to the four-fermion 



interaction) to the muon decay amplitude. According to Sirlin p3 

a ,„ 7 -As 2 



° = 4^ (6 + ^f- lnc2 + 4Aw) ' (K8) 
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where 

2 

A w = A(m w ) = j—^ + ln4vr - 7 - \n(m 2 w /fi 2 ) (E.9) 
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Appendix F 
The radiators R,A q and Rvq- 



For decays to light quarks q = u,d, s, we neglect the quark masses and take into account 
the gluon exchanges in the final state up to terms ~ - pE5|] , and also one-photon 



exchange in the final state and the interference of the photon and the gluon exchanges [[19[ 
These corrections are slightly different for the vector and the axial channels. 
For decays to quarks we have 

T q = T(Z - qq) = l2[g 2 Aq R Aq + g 2 Vq R Vq ]T (F.l) 

where the factors Ra,v are responsible for the interaction in the final state (in our previous 
papers, we used letter G instead of R) . The results presented in (45| - (49| are 



Rv q = 1 + — + -Q \--- A Q -- + 

71 4 y 7T A H 7X 7X 



l n ,a. SN9 ,cL,, . „„m 2 d. 



+ [1.409 + (0.065 + 0.0151nt)^](^) 2 -12.77R) 3 + 12-f-^<U (F.2) 



t ix n m z 7X 



^ = R Vq - (2T 3q )[I 2 (t)(^Y + h(t)( 

7T 7T ' 

^d s m 2 1 m 2 a 2 2 

" 12 ^l T*™ " " 10 W ( V } *™ (R3) 



where m g is the running quark mass (see below), 



r OL x 



8.7(— )+45.15(— ) 2 , (F.4) 

7T 7T 



C = 1 + 3.67(2!) + (11.29 - lnt)(^) 2 , (F.5) 

7X 7X 



, 2 ( iH -3.083 -ln i + ^? + ^ , (F.7) 



J 3 (i) = -15.988 - 3.722 hit + 1.9171n 2 1 , (F.8) 
t = m 2 /m| . 

Terms of the order of (^-) 3 caused by the diagrams with three gluons in intermediate state 



were calculated in |79[. For Ry q they are numerically very small, ~ 10 5 ; for this reason, we 
dropped them from formula ( F.2|) . 
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For the Z — » bb decay, the 6-quark mass is not negligible; it reduces Tb by about 1 MeV 
(~ 0.5%). The gluon corrections result in a replacement of the pole mass m& ~ 4.7GeV by 
the running mass, the virtuality being mz ■ — > rhb(mz)- We express rhb{mz) in terms 
of nib, « s (^z) and a s {mb) using standard two- loop equations in the MS scheme (see |50|). 

For the Z — > cc decay, the running mass rh c (mz) is of the order of 0.5 GeV and the corre- 
sponding contribution to r c is of the order of 0.05 MeV. We have included this infinitesimal 
term in the LEPTOP code, since it is taken into account in other codes (see, for example, 
0). 

We need to remark in connection with T c that the term ^{t), given by equation ( [F.7| ), 
contains interference terms ~ (a s /ir) 2 . These terms are related to three types of final states: 
one quark pair, a quark pair and a gluon, two quark pairs. This last contribution comes to 
about 5% of I2 and is infinitesimally small at the currently achievable experimental accuracy. 
Nevertheless, in principle these terms require special consideration, especially if these quark 
pairs are of different flavors, for example, bbcc. Such mixed quark pairs must be discussed 
separately. 

Note that a s stands for the strong interaction constant in the MS subtraction scheme, 
with /x 2 = m 2 z . 



69 



Appendix G. 

Derivation of formulas for the asymmetries. 



Asymmetry in the processes e + e~ — > Z — > ff is calculated with the masses of e and / 
neglected in comparison with the Z boson mass. (Mass corrections for f = b are of the order 
of 2 • 10~ 3 and will be taken into account below). The amplitude (20) of the interactions of 
the Z boson with massless fermions // can be conveniently rewritten in the form 

M(Z -> ff) = \fZ a W L 3« + 9 f RJ«] , (G.l) 



where 



jt R = ^L,Rla^ 
^L,R = -(1±75)^ , 

9l,r = 9vf±9Af ■ 

The chirality is a conserved quantum number for massless fermions (anomalies do not yet 
manifest themselves in the approximations we deal with here) and coincides with a fermion's 
helicity up to the sign. 

Therefore, the pairs t^t\ and e R e R do not transform into the Z boson at all, and the 
pairs e^e^ and e R e~j^ create a Z boson with the polarization ±1, respectively (along the 
positron beam). The scattering amplitudes thus have the form 

T (eZ,i?e + -> f L ,Rf) = #2,/?£l,/? t o(1 + cosfl) 
T(e- LiR e + - f R , L f) = g e L>RgR>L T (l - costf) , (G.2) 

where coefficient T is nonimportant at the moment (it can be reconstructed from ( |G.1| )), 
and 8 is the angle between the momenta of e~ and /. The sign in front of cos 9 is chosen for 
the helicity to be conserved in forward and backward scattering. 

Once the form of the amplitude ( |G.2| ) is known, all asymmetries are immediately found. 

(a) Left-right asymmetry Alr is defined as the ratio 

, _ Ql- or 

A LR — : , 

0"L + OR 

where a L<R = a(e LtR e + -> //). Hence 

A (9l) 2 ~ (9r) 2 _ A e t n o\ 

AlR ~ ( e\2 i f e\2 = A V^- 6 ) 

Wl) 2 + [9r) 

(b) Longitudinal polarization P T (cos#) is defined as the ratio of the difference to the sum 
of differential cross sections, = ^f(ee — > Tr^t): 

( da\ / da \ 

PA<**0)= )Z\\ )Z\ L , (G.4) 



L,R 
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where 



,da 

da, 
{ dd' 



2m| 



2m? z 



To\\9 R Y[{g R Y{i + ^ey + {giY{i- cos 



ro|W[^ir(l + COB^ +^(1-008 



(G.5) 



Substituting ( G.5|) into the definition ( |G.4| ), we obtain 

A T (1 + cos 2 9) + 2A e cos 9 



P-rfCOS I 



'1 + cos 2 6 1 ) + 2A e A T cos 9 



(G.6) 



where A e and A r are defined according to ( |G.3| ). The longitudinal polarization P T , averaged 
over directions of r-leptons, is defined as the following ratio of the total cross sections <tl r = 
J\d cos 9(%) LtR 



°r ~ °l I-i dcos9[A T (l + cos 2 9) + 2A e cosi 



°-R + <r T L 
A T 



j\dcos9[l + cos 2 9 + 2A e A T cos 9} 



(G.7) 



(c) Forward-backward asymmetry Ap B is calculated more simply in terms of gA,v- The 
squared matrix element of the process ee — > ff is proportional to 



\M\ 2 ~ {[g 2 Ae + g 2 Ve ] [[g% + g 2 Vf ](l + v 2 cos 2 



+ [9 2 Af-9 2 Af\(l 



+ ^(gVegAegVfgAf)v COS 9 



(G.8) 



4m / 



where 9 is the scattering angle and v — 1 — -^J- is the velocity of fermion /. This immediately 
implies that 



A f 



FB 



ZgAfgvfv 



9 2 AfV 2 + g 2 vf 3 -^ 



(G.9) 



The mass m/ is negligible in all channels with the exception of / = b, where the nonzero 
mass produces effects of the order of 2 • 1CT 3 . Gluon corrections in the final state (see 
Appendix F) replace the pole mass ~ 4.7 GeV in equations ( |G.8| )-( |G~9| ) by the running 
mass at the mz scale: — > rhb(mz)- 

Note that starting from gluon corrections of the order of (— ) 2 , it is impossible to un- 
ambiguously separate different quark channels, since additional pairs of 'alien' quarks are 
created in this order. We do not consider corrections (^) 2 in asymmetries. In our approx- 
imation the ratio gvf/gAf is n °t renormalized by the gluonic interaction in the final state. 
Therefore, the expected accuracy of formula ( |G.9| ) is (^) 2 ~ 2 • 10 -3 , which is by an order 
of magnitude better than the experimental accuracy. 
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Appendix H. 
Corrections proportional to m^. 

This appendix gives a simple mnemonic recipe for the derivation of corrections propor- 
tional to m 2 . A rigorous derivation requires careful regularization of Feynman integrals. 

The terms proportional to m\ contribute to radiative corrections to bare masses (squared) 
of the W- and Z-bosons, but not to the corrections to the bare coupling constants. This 
follows from dimensional arguments. Indeed, the dimension of self-energy E for the boson, 
equals m 2 ; therefore, the terms ~ m 2 remain in S(g 2 ) in the limit q 2 — > 0. On the other hand, 
the corrections to coupling constants are proportional to dH/dq 2 and do not contain terms 
~ m 2 . Therefore, it is easy to evaluate the contribution of the t-quark to the parameter 
P — ( a z I '&w)( m w l m z) i n the approximation ~ am 2 (the Veltman approximation 
neglecting the terms ~ a: 

_ ®zo - £ w (m^) _ ^ S z (0) T, w (0) _ 
a wo m| - S z (m|) ~~ m| m 2 w 

= i + n z (o)-n w (o) . (H.i) 

The evaluation of the difference ITz(O) — Iljy(O) is elementary: 

5V(0) 



nv(o) 



m 



w 



3a w ( f°° p 2 dp 2 f°° 2 2 f°° dp 



Snmly \Jo p 2 + m 



rl dp2 -^r?r^) (a2) 



(As we have neglected the mass of the 6-quark, the propagator of the 6-quark compensates 
the factor p 2 in the numerator.) 



IWO) - ^ 



rrir 



3a z fl [°° , 2 1 [°° P 4 dp 2 
\2Jo P + 2Jo 



]7im 2 z \2 Jo LJ 2 Jo (p 2 + m 2 ) 2 

lo p z + mf 

Taking into account that in one-loop approximation we can set in ( |H.2| ) and 

(H.4) 



aw otz 



m z 

we see that quadratic and logarithmic divergences cancel out and that finally 

3ctz m 2 „ 3az /TT _ % 

p ~ 1 + Ap t = 1 + — Z-L = l + —^t, (H.5) 

107T m z 107T 
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where t = m|/m| and we assume that t ^> 1. 

Let us express the leading in t corrections to the main quantities m^/m^, and gy l in 
terms of Ap t . 

We define 

4 = a w /a z ,4 = 1 - 4 ■ (H.7) 

Then 
and hence 

44 * s 2 c 2 /(l + Ap t ) . (H.9) 
Solving the last equation, we obtain 

4^ (l + ^A Pt ) , (H.10) 



c 2 



and therefore, 



and 



Likewise, 



<-* 2 U--T— ^ > (H-ll) 



* 4(1 + A*) ^ c 2 f 1 + , (H.12) 

m z \ c z — s z J 

^~ C+ J5. £ t mi3) 



1 - 4s' 



: (1 -4s 2 ) + 
47r(c 2 — s 2 



4c 2 s 2 
c 2 — s 2 



Apt 



(H.15) 



The corrections proportional to m 2 were first pointed out by Veltman p9|, who em- 



phasized the appearance of such corrections for a large difference m 2 — m 2 which violates 
the isotopic symmetry. In this review, the coefficients in front of the factors t in equations 
( |H.13| )-( |H.15 ) are used as coefficients for normalized radiative corrections Vi. 
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Appendix I 

Explicit form of the functions T\(t) and Hi(h). 



The equations for Ti(t) and Hi(h) are [pffi, [63|: 



i = m 

m , . .2 8 n. , 4 32 o 

T m (t) = s 2 ) hit - - + — s 2 + 

w v 3 9 ; 3 9 



2c 2 - s 2 2 lc 2 - s 2 r 2 16 2 2 32 2 lri/ 

3 — t 2 + - — t + s 2 - -t s 2 t ]F t (t 

3 c 4 3 c 2 L 3 9 3 9 J V 

rr /, s h . . c 2 h ^ h s 2 , 8 2 

# m (/0 = -T 7 In h + 2 In - - — -/i - -s 2 + 

h — 1 h — & c A 18<r 3 



> 2 4/i 4 Nri/lN 

+ <¥-T + 3 )F " (ft) - 

+ (1.1203 - 2.595s 2 ) ; 
where 5s 2 = 0.23110 — s 2 (note the sign!). 

i = A 

T A (t) = I - S -s* + - 1 - 2ff '"> + 
AW 3 9 27 At - 1 

,32 4 8 2 l Nr 4^ /x 2,, „.l-2tF t (t). 

+ (— s 4 - -s 2 - -) -tF t (t) - - 1 + 2t) ^ 

v 9 3 2 ;[ 3 W 3 V ; 4t - 1 J 

, n c 2 . /i 8/i . , 

-—h+ [0.7752 + 1.075s 2 ] . 
18 

i = R 

T R {t) = 2 -\nt + ^- 2 -{l + llt)F t {t) ■ 

Ah c 2 h 

H R (h) = 1 Try- r In -7 + 

Ry 1 3 18 1 - c 2 /h c 2 
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+ <!-§k + i*')ft W + T A s b.* + 

+ (1.3590 + 0.5L5s 2 ) 
i = v 
T v (t) = T A (t) 

H u {h) = H A {h) (1.4) 

The functions F t and are the limiting cases of the function F(s, mi, m-i), described in 
Appendix B. The explicit formulas for F t (t) and Fh(h) are given in Appendix J (equation 
(J. 3)) and Appendix L (equation (L.ll)), respectively. 
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Appendix J. 

The contribution of heavy fermions to the self-energy 

of the vector bosons. 

Let us give the expressions for the contribution of third-generation quarks (t, b) to the po- 
larization operators (self-energies) of the vector bosons. We use the following notations: t = 
m 2 /m 2 z , b = m 2 b /m 2 z , (b < 1), h = m 2 H /m 2 z , IT 7 (g 2 ) = £ 77 (g 2 )/m|, Ii lZ (q 2 ) = £ 7Z (g 2 )/m|, 
U w (q 2 ) = Z w (q 2 )/m 2 w . 

The dimensional regularization yields the terms 

2 m 2 
A, = — — - 7 + ln4vr-ln^- , (J.l) 



where i = t,b,W, Z..., D is the variable dimension of spacetime, (4 — D = 2e, e —>■ 0), 
7 = -F(l) = 0.577. . . (we follow |0|, p. 53). 

We begin with an auxiliary function F t (t), obtained as a limiting case of the function 
F(s,mi,m 2 ) (see Appendix B and [|(J, p 54; ||8l| , p 88), 



F t (t)=F(s = m 2 z ,m t ,m t ) = F(l,t,t) , (J.2) 



and get, using f8~0fl , 

| 2[1 - y/At- larcsin^] , At > 1 , 
Ft{t) = { 2[1 - VT^Atln i±$p , 4t < 1 . (J - 3) 

The asymptotics of F t are 

F t ~lni + 2 t->0 , 

F t ~ l/6t + l/60t 2 t -> oo . (J.4) 

Differentiation gives 

2 dF d „ 1 - 2£F t 

F/ = m 2 — r = -t— F t = . (J.5) 

4 z dm% dt At -I K J 

In this Appendix, IT stands for the contribution of the doublet (t, b) to the corresponding 
polarization operator: 

n 7 (o) = o . (j.6) 



n 7 (m|) = - 

7T 



Q 2 ^A t + (l + 2t)F t (t)- 1 -^+Q 2 b ^A b + ^ + \nb^ . (J.7) 
n 7Z (o) = , (J.8) 

n 7 z(m|) = ^-{{9l-s 2 Q 2 )[A t + {l + 2t)F t {t)- 1 ^- 

Q b , _2/->2\ ( A , ^ 



- If + s z Qi)[A b + - + \nb)\ . (J.9) 
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iMo) 



a 



Atxs 2 c 2 



3 . 3 
2* A * + 4* 



(J.IO) 



a 



4tts 2 



3t \ 5 * 

2^2; ' + 3 ~ ^ 



2c 4 ' 



2c 2 + 2c 6 



In II - 



(J.ll) 



nz(m|) 



+ 



as 

7TC 2 

a 

2^? 2 L 
a 

871S 2 C 2 



Q 2 t (A* + (1 + 2t)F t (f) - i) + Q 2 (A, + I + In ft) 
Qt (A t + (1 + 2t)F t (t) - - g b (A, + I + In 6) 



(2-3t)A t + (l-t)F t (t) + - + lnt 



n z (o) 



a 



Ans 2 c 2 V2 



-tA, 



+ 



(J.12) 
(J.13) 



E'(m|) 



+ 



+ 



as 

7TC 2 

a 

27? 

a; 



Qt (a 4 + F t - - + (i + 2*)f;) + q 2 (a 6 + - + In ft) 

Q t (A t + F t - 1 + (1 + 2t)F/) - Q 6 (a 6 + + lnfc) 



87TS 2 C 2 



2A t + F t + - + ln* + (1-0^-1 



as 2 



TVS" 



Qt A, + 



3(4t - 1) 



2t 1 
+ - 



2t + 4t z 



+ 



a 



2ttc 2 
a 

8tts 2 c 2 



1-4* 
l-2t + At 2 



2 + t 

2A ' + lnt + 3(¥3T) + 



1-4* 
2t 2 + 2t - 1 



4t- 1 



+ 

(J.14) 



nz(m|) 
a 



3a 



z) -= [2tF t (t) - (1 + 2t)F/(t)][^g 2 - ^- 2 Q t 



167TS 2 C 2 87TS 2 C 2 



a 



a 



as 2 



167TS 2 C 2 



(2-3F t '(t)) + — Q 6 + — Q 



2ttc 2 



7TC 



2 

2^6 



(J.15) 



Substituting the expressions for the polarization operators into the formulas for physical 
observables, we verify the cancelation of the terms ~ A$, and also the terms proportional 
to In b, since the limit m h — > does not produce divergences. It is convenient to get rid of 
the terms ~ A b and ~ In 6 already in the expression for the polarization operators using 
equation 

A 6 + ln(6) = A t + ln(t) (J. 16) 
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and then making sure that the terms ~ A t are indeed eliminated. 

Our definition of the wavefunction of the Z-boson differs in sign from that assumed in 
80] ; hence the quantity U^ z we use also differs in sign from the expression given in |80|| . With 
our definition, the interaction of Z bosons with Weyl fermions is —ief^y^fiTs — Qs 2 )Z fl , and 
that of photons is —ieQf'j^fA^. The latter vertex coincides with the one given in equations 
(8)-(9) of ||80|| , while the former differs in sign. 

Let us look at the formulas for physical observables. 

The quantity T m (t) (|I.1|) is defined as the following combination of polarization operators: 



t + T m {t) 



IQtcs 4 r c 2 



— {-[U z (m 2 z )-U w (m 2 w )] + 



3a 

+ U w (m 2 w )-U w (0)-U^m 2 z )}, 



(J.17) 



since I1 7 ^(0) = for fermion loops. 

Using equations (|J.7|) , ( |J.10|) -( |JTT2"[ ) we obtain 



t + T m (t) 



t+[--~s' 



+ 



+ 



(c 2 - s 2 ) 



4 

c 6 c 2 , 



32 
— + 
9 

In II - 



c 2 -s 2 



c 



3c 2 t + 
2 ^2(c 2 - 



t 



3c 4 



-t 2 + 



16 

IT' 



2 -t-fsH)F t (t) 



(J.18) 



where we have taken into account that Il 1 {m 2 z ) cancels out with the sum of terms propor- 
tional to Q 2 and Q 2 in Uz(M z ). The terms proportional to t arise from II w (0) and T1 W (M Z ); 
those proportional to t 2 arise from Iiw{ m w) \ the terms proportional to Int and the constants 
arise from n^y(m^ / ) and U z (m 2 z ); In |1 — y| corresponds to the threshold (tb) in liwi.'^w)^ 
and finally, the only source of terms proportional to F t (t) is U z (m z ). 

The infinities in the contribution of the doublet (t, b) to the observables must cancel each 
other since the introduction of an additional fermion family into the electroweak theory does 
not violate its renormalizability. 

Substituting the terms proportional to A t in ( |J.17p and taking into account (|J.16|) , we 
obtain zero: 



a 



8ns 2 c 2 



(2 - 3t) - 



a 



2irc 2 



+ 



a 



Airs 2 



3tA 
2c 2 ) 



Airs 2 
+ 



1 



Atts 2 c 2 2 



3t 

> 

V 



+ 



The expression for T^(t) is 
t + T A (t) = 



16tts 2 c 2 
3a 



U z (m 2 2 



- E' z (m|) - n 



(0)] 



Using (pTT0|) , (|JT2l) and ( PH) , we have 

2 



t + T A (t) = t + 



8 2 , 16 4 
-s H s 

9 27 



F! + 



32 



s 4 - 



x 



2t) 



(J.19) 



(J. 20) 



(J.21) 
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The terms proportional to At obviously cancel out. IIvk(O) gives a contribution propor- 
tional to t, while all other terms arise from the difference Uz{m 2 z ) — T! z {m 2 z ). 
Finally, we look at T R (t): 



t + T R (t) 



1Qtcc 2 s 2 
3a 



(c 2 - s 2 : 



cs 



■n Z7 (m|) + n 7 (m|)-n z (m|) + n w (o)l . (3.22) 



The terms proportional to Q 2 and Q 2 cancel out, so that only 11^(0) remain, as well 
as terms that do not contain Q 2 b coming from ITz 7 and liz- Substituting ( |J.10| ), (|J.9|) and 
QJ.12D , we have 

t + T R (t) = t + ^ + ~lnt-~(l + llt)F(t) . (J.23) 

Here t comes from IIp^(O), the term proportional to F t {t) - from IT^, and | + §hit - 
from I! z and n^ 7 . 

Cancelation of infinities in flj.22|) follows from 



A t 



I 2 2^ 
C — S 



a 



Aires 



cs 



a 



(2 - 3t)a 3ta 



2nc 2 



8nc 2 s 2 



8tts 2 c 2 



. 



(J.24) 
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Appendix K. 



The contribution of light fermions to the self-energy of 

the vector bosons. 



The contribution of the doublet of light fermions to the polarization operators is readily 
obtained using the formulas of the preceding Appendix. 

To achieve this, Az,w must be substituted into the formulas of Appendix J instead 
of A q + ln(m q /mz t w) 2 ■ The physical reason for the absence of terms proportional to the 
logarithm of the mass of light quarks and leptons is the infrared stability of the quantities 
that are analyzed in this Appendix. In the formulas (|K.4j ) and ( |K.7|) we use the equality 
Qu — Qd = 1- The subscripts u and d stand for the upper and lower components of the 
doublet: 

n 7 (o) = o . (k.i) 

n 7 (m|) = ^(Ql + Ql)(A z + |) . (K.2) 
IV(0) = . (K.3) 
n 7 z(m|) = ^{A z + (Q5 + Q> 2 ] . (K.4) 

oCSTT O 4 

nv(o) = o . (k.5) 

^w(m 2 w ) = ^- 2 (A w + 5 -); . (K.6) 

n ^ (m|) = S^^ z s ^ Ql + Q2 * )] ■ (K - 7) 

n z (0) = . (K.8) 

^^ = ^(^ + l)l\-Y + ^Ql + Q 2 d )] ■ (k.9) 

Formulas ( |K.lj )- [KI^ ) must be used for three lepton doublets (y e ,e), {v^, fJ>) and (y T , r) 
with N c = 1 and two quarks doublets (u, d) and (c, s) with N c = 3. 

Substituting equations ([K.I]) - ( |K.9| ) into expressions for physical observables in terms 
of polarization operators (|L.16|) , (|L.20| ) and ( L.24[) , we arrive at the contributions to the 
constants Cj owing to the self-energies. 
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Appendix L. 

The contribution of the vector and scalar bosons to the 
self-energy of the vector bosons. 



This Appendix gives formulas for boson contributions to polarization operators that we 
reproduced from |[C], pp 53, 54. (There is a misprint in |8(|: the term proportional to A\y 
in the expressions for Uw(q 2 ) must be multiplied by 1/3). 

The polarization operators in ]50| depend on cw and sw via coupling constants and de- 
pend dynamically on the ratio mw/n^z, which arises from Feynman integrals. We substitute 
everywhere c for cw (and mw/mz) an d s for and sw In the framework of the one- loop 
approximation, this substitution is justified. After this substitution, we find expressions for 
physical observables; ultraviolet divergences of polarization operators cancel out in these 
expressions. 

In the following formulas Ilj denotes only boson contributions to the corresponding po- 
larization operator (all calculations were performed in the 'tHooft-Feynman gauge): 



n 7 (m|) 



a 



^^3A^ + 2(3 + 4c 2A 
Air 



1 - V4c 2 - 1 



arcsm 



a: 
Air 



(3A W + 1.53) . 



n 7 (o) = o . 



n 7 z(o) 



a 



Aires 



(2c 2 A w ) 



ca 



A 



w 



2c 



(L.l) 

(L.2) 
(L.3) 



n 7Z (m|) 



a 



{(5c 2 + i)A 



Aires 

A 13 2 
+ 2 - + — c 2 
V 6 3 

-a /30c 2 + 1 
Air 1 6cs 



4c 4 



A 



1 

2c 



1 — V4c 2 — 1 arcsin [ — 



3.76 



(L.4) 



iMo) 



a 



Airs 2 

„4 



' s 2 \ 3 c 2 

~ 2 ~ l ) Aw + A{l-c 2 /h) ln l 



tr 1 39 s 2 



a 



Airs 2 



8c 2 12 

i ~ l ) Aw+ 



h 

8c~ 2 



In — - 



h 



A(l-c 2 /h) h 8c 



+ 0.85 



(L.5) 
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+ 



+ 



+ 



+ 



4ns 2 



1 + 17-2 3 



c 4 



2c 4 J 



cr 1 3 \ s 4 1 
- - + 21 -- + - 



2 1 s 4 113 



3c 2 12c 4 



F,(Vc 2 ) - 1 - — + 



6c 2 
3h 



Qc 2 4(c 2 -h) 



18 
\n(h/c 2 ) 



a 



4ns 2 



25 



h 



___ Aw _ 1 . 76+ M__ 



3/iln(/i/c 2 ; 
4(c 2 - fe) 



h 2 \ ^ , 9x /i 
12? Wc) -6^ 



nz(m|) 



+ 



x 



+ 



7s 2 - 



25 



4tts 2 [ \ 

1 + 8c 2 



7s 2 \ 73 2 

6? Aw + 36?"9 



2 + — r^-(c 2 - s 2 ) 2 -yc 2 (l + 2c 2 )^ 



6c 2 

1 — V 4c 2 — 1 arcsin ( — 
V2c 



h 



c 6c 2 4c 2 (1 - /i) 



~2 

a 

47TS 2 
/ 1 



+ 



3h 



\nh + 



h 



+ 



13 
12c 2 



ln(c 2 ) 



/i 2 \ 



7s 2 
/i 



— + --] A 

6 6 c 2 / 147 



; 2 3c 2 12c 2 
0.58- A + 



3h 



Qc 2 4c 2 (l-/i) 



In /i 



3c 2 + l£) ^ 



M = 



+ 
+ 

x 

+ 

+ 



ds U=m2 z 
— ' 3 - — + — ) A w 



4n 



a 



6s 2 6c 2 / 
40c 4 



C 2 - 5 2 ) 2 
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12c 2 + (8c 2 + l)(c 2 - s 2 ) 2 - 40c 4 (l + 2c 2 ) 
4c 2 / 1 

-1 + 



1 — V 4c 2 — 1 arcsin 



2c 



= arcsin , 
V4c 2 - 1 V2c 

1 - (h - l) 2 ] F h (h) + [ll - 2h + (1 - h) 2 } F' h {h) 

1 • // , , 1 h 1 2 

1 ~ 2(fa- 1) r ~ + " [ + (r ' s ) i! " 
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4vrs 2 c 2 V 3 12 hK J \6 12 " v ' 



+ 



h 



Inh- 1.67 



12(1 - h) 

The functions Fh(h) and F^h) are defined as 

F /t (/i) = F(s, m z , m H ) \ s = m %= F(l, 1, h) , 
dF(s,m z ,m H ) 



F' h {h) = s- 



ds 



Using |pl] , pp 88, we obtain 



(L.8) 

(L.9) 
(L.10) 



F h (h) 



(L.ll) 



If /i — > oo, 



If h -> 0, 



1 + I ) In /i — /ia/ - — 1 arctan \ 1 , h < 4. 

V h — 1 2 / V/i V/i 



r, 1 1 A 4 \ 5 59 

F ^^2^-^( 1 + ^) ln/i+ 3^ + T2^ 



» 1 -nVh + (1 - -ln/i)/i . 



(L.12) 
(L.13) 



Finally, for F^(h) we have 

h-1 



-1 



2 ln ' + (3 "W/^ ln vV + Vl 1 ' /,:> 1 



(L.14) 



If /i — > oo, 



h-1, , , ,, /i /4-/i 

-1 H In /i + (3 — h)\ arctan \ , h < 4. 

2 v '\A-h V /i 



(L.15) 



All infinities in formulas ( p^ . 1| ) — (|L78|) are collected into the factors A^ by replacing the 
factors A z and A H using the equation 

A; = Aj + ln(m 2 /m 2 ) . 



The function H m (h) is 
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+ H w (m 2 w )-H w (0) 

- n 7 (m|) - 2-H^ z (0)) - divH m , (L.16) 

where divH m denotes the sum of terms, proportional to Apy, in polarization operators in 
( |L.16|) . Substituting the finite parts of the formulas for polarization operators from this 
Appendix, we obtain 



H m (h) 



1U /I s 

In 



^-z{m 2 z ) - H w (m 2 w ) 



H m {h) 



h 



h-1 



- (c 2 -s 2 ) 
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9c 4 



c 2 h 


111 2 

C 2 


s 2 


h-c 2 


18c 2 


Ah 






~ 9? 





x 9 

0.50 



Ah A 
~9~ + 3 



(L.17) 



where the term proportional to Fh(h) arises from Hz, and the term proportional to F^h/c 2 ) 
arises from Hw(m w ). The term proportional to Inh originates from Hz, while \n(h/c 2 ) 
arises both from H w (m 2 v ) and from 11^(0). The term proportional to h is contained in H z , 
^-w{ m w) an d nw(0) and finally, all four polarization operators make contribution to the 
constant. 

Collecting the terms proportional to Ajy in the polarization operators in ( |L.16j ), we get 



divHr, 



+ 



167TS Z 

3a 
a 
Ans 2 



A 



w 



a 



25 7 s 



25 



Ans 2 s 2 



7s z - — + -— + 



6 



19\ 3a sea 
— + — + 2— — 
by Air c s 2tt 



6 c 2 6 

16vrs 4 . 
A 



a 



3a 



irs' 



(L.18) 



Note that the divergent term in D (see (|E.8j ) and ( |M.26|) ) exactly compensates for the 
divergence in ( L.18|) , which justifies the subtraction of infinity in ( PL .16 ). 

The function HA{h) is expressed in terms of polarization operators as follows: 



H A (h) 



16irs 2 c 2 



3a 



H z (m%) - E' z (m|) - H w (0)\ - divH A 
Substituting the finite parts of the polarization operators, we obtain 

8h 



(L.19) 



H A (h) 



he 2 , h 
— In — 



h — c 
(A 



9(h 



-Hh) + ^-h + h 2 )F h (h) 



Ah h 2 \ h 

+ — F'(h) + 0.78 

3 9 9 } hK J 18 



(L.20) 



where \n(h/c 2 ) stems from ITvv and In h stems from Hz- Fh(h) arises both from Hz and from 
TI Z , while the only source of F' h {K) is TI Z . The term linear in h is contained in Hz and Hw, 
while all three polarization operators contribute to the constant. 
Adding up the divergent terms, we have 



16 

divH A = —c 2 s 2 A w 



16ns 2 c 2 . a 
3a 7r 



(L.21) 
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Note that the divergent part of D A (see ( [jVLlgj) and (|M14|) ) is: 

16c 2 



divD, 



-A 



167TS 2 C 2 . Q! 

-A w — 



3 3a 7ts 

Vertex parts also contain ultraviolet divergences (see Appendix M, eqs.( M.iT ), ( |M.15| )): 



(L.22) 



cUvFa 



16 



-c s 



16 
~3~ 



2„2 - C 2 )A 



167rs 2 c 2 c 2 a 
w — — ^— — A w 



3a 



S 2 7T 



(L.23) 



The sum of terms (|L.21|) , ( [L.22|) and (|L.23| ) equals zero. 
Finally, the expression for H R (h) is 



H R (h) 



16tt 
3a 



2 2 
C S 



(c 2 



C'.S 



-n Z7 (m|) + n 7 (m|) - n z (m|) + n w (o) 



+ 2-n 7Z (0) ) - divH R (h) . 



Collecting the finite parts of the polarization operators, we find 
H R{ h) = -^ + r ^ In ^ + (| - ^ + I/, 2 ) ^) + ^ In A + 0.03 . 



(L.24) 



(L.25) 



The term proportional to Fh(h) stems from Hz, just as In h does. generates the term 
~ ln(/i/c 2 ). The term linear in h is contained both in IT Z and in Tlw, and all polarization 
operators with the exception of IT 7Z (0) contribute to the constant. 

Adding up the divergent parts of the polarization operators in ( L.24Q , we obtain 



divH R (h) = A 



IV 



^(c 2 - s 2 )(l + 30c 2 ) + 4cV + ^c 2 ( 7s 2 - ^ 
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(L.26) 
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4 2 ^ 2 i 16 2 2 



-s H c s 
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Aw(8c 2 - ^c 4 ) 



Taking into account the divergent term in D, yielding (see eqs. 

16 



divD 



R 



-c 2 A 



IV 



), prm ( To) ) 

(L.27) 



and in F R (see Appendix M, equations (M.ll) and (M.22)), 

16 



divFf, 



-c 4 - =-c 2 )A w 



3 3 

we confirm that divergences cancel out in expressions for R. 



(L.28) 
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Appendix M 

The vertex parts of FAf and Fyf and the constants Ci. 



This Appendix collects the vertex functions that form a part of one-loop electroweak 
corrections to Z — > up, Z — > Z + / _ , Z — > mm, cc, c?(i and ss decays. In the case of the Z — > bb 
decay, a t-quark can propagate in the loop, so vertex corrections are not reducible to numbers 
but are functions of m t (see Appendix N). 

The finite parts of vertex functions are given in [3(J, pp 29, 30. The corresponding ex- 
pressions depend on cw{sw) and mw/mz- In the framework of the one-loop approximation, 
we replace cw and mw/mz with c, and sw with s. For this reason, while vertex functions 
in [BD|| depend on m t , itlh and the new physics, ours are numbers (see also Appendix L). 

This Appendix also gives the infinite parts absent from Q and required for testing 
whether the infinities in physical observables do cancel out. 

We begin with the Z —>■ vv decay: 



F, = F\ 



Vv — 

a 1 
Att Acs 



1 Au 



-A 2 (m 2 z ,m z ) + 



2s 2 



4c 2 s 2 *\ *> »/ 2S 2 

For the decay to a pair of charged leptons or quarks we have 

v f (v 2 f + 3oJ)A 2 (m|, m z ) + F f L 



1 3c 2 
-A 2 (m 2 z ,m w ) + — A 3 (m 2 z ,m w ) 
s z 



F - a 



(M.l) 



(M.2) 



P - a 

FAf -A^ 



a f (3v 2 f + a))A 2 (m 2 Z: m z ) + F[ 



(M.3) 



where a\ 
Q l = -1, Q 



a d = -1/Asc, a u = 1/Asc, v f = (T/ - 2Q^s 2 )/2sc (Tj = T 3 d = -1/2, T 3 " = 1/2, 
d - -1/3, Q u = 2/3). The functions F[ are 



Ft 



Ft 



1 3c 

— A 2 (m 2 z ,m w ) - —A 3 (m 2 z ,m w ) 



-s 2 3c 
—A 2 (m 2 z , m w ) + — A 3 (m|, m w ) , 



]s 3 c 



1 - ± s 3c 

A 2 (m|, m w ) + — A 3 (m|, m w ) 



(MA) 

(M.5) 
(M.6) 



For calculating Fyf and Fa/ we need to determine the values of three constants: A 2 (m 2 z , mw) 
A 2 (m z ,m z ) and A 3 (m z ,m w ): 



A 2 (m 2 z ,m w ) 



--- 2c 2 -(2c 2 + 3) ln(c 2 ) 



+ 2(1 + c 



2\2 



In c 2 In f — ^_ — 



S P -3 



(M.7) 
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where we have used mw/mz = c; Sp(x) is the Spence function: 

\ (It 7T 2 

Sp(x) = -Jj ln(l - art) , Sft-l) = -- . (M.8) 
o 

Using (|M7?|) and (pO|) , we find 

A 2 (m|,m z ) = -- - 2 - 85-p(-l) . (M.9) 



Finally, 



A3(m z ,miy) = c H — (2c + l)v4cr — 1 arctan ■ 



6 3 3 V ' V4c 2 - 1 



3— '-r^Tl^T) ■ (M - 10) 



— -c 2 (c 2 + 2) I arctan 



The expressions for divergent parts of the vertex functions, describing the coupling of 
the Z-boson to the leptons, are 

a c 

divF v = n^-w ,divF V i = divF A i = —divF u (M.ll) 

87T s A 

We switch to the calculation of the constants C{. We begin with definitions. According to 

V v {t, h)=t + T„(t) + H v {h) +L V + D U + F U , (M.12) 

The value of L v represents the contribution of leptons and light quarks to the polarization 
operators of the vector bosons and can be easily obtained from the formulas of Appendix K 
for polarization operators: 

L v = 4 - 8s 2 + ^s 4 . (M.13) 

D u originates from the box and vertex electroweak corrections to the /i-decay The 
expression for D see in Appendix E, equation (E.8). For D u we have 



16tts 2 c 2 n , 
D„= ^—(D -A w ) . M.14 

6a 7TS 



Finally, 



h = ™g* F . . (M.15) 
6a 

Comparing ( [IVl . 1 2| ) and (|58"D, we arrive at the expressions for C v whose ingredients are now 
all determined: 

C U = L U + D V + F v (M.16) 

Let us switch to V A : 

V A (t, h) = t + T A (t) + H A (h) + L A + D A + F A . (M.17) 
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The expressions for La and Da are already there: 

L A = L V , D A = D V ; (M.18) 
the formula for the vertex function is 

h _ j*gt FM . (M . 19) 

Finally, 

Ca = L a + D a + F a (M.20) 

We now move to V R : 

V R (t, h)=t + T R (t) + H R (h) +L R + D R + F R , (M.21) 

where D R = D v and Lr = since n^(m^) is absent from the ratio gv/dA- 
For F R we have 

= 16?r(c2 " 82)08 [-JVi + (1 - 4 S 2 )F^] , (M.22) 
6a 

and the formula for Cr is 

C R = L fl + D R + F R . (M.23) 
We end this Appendix with formulas for C m : 

C m = L m + D m , (M.24) 

L m = 4(c 2 -s 2 )lnc 2 , (M.25) 

D m = - l -^(D-^ w ) . (M.26) 
3a V ns 2 J 



88 



Appendix N. 

The functions 4>(t) and S(f>(t) in the Z — > bb decay. 

For the function <p(t) we use the expansion from R57|| : 

3 -2s 2 f , , t 



2s 2 c 2 



\t + c 2 \2.mn— -6.716+ 



+ -(8.368c 2 In 4r - 3.408c 2 ) + 4(9.126c 4 In 4r + 2.26c 4 ) + 
t c 2 t 2 c 2 

+ ^(4.043c 6 In + 7.41c 6 ) + ...} , 
and for 5<f)(t) we use the leading approximation calculated in ]5£j and 

1 



(N.l) 



... , s 3 — 2s 2 f 7T 2 cL(m t ). 



(-)tVf(^) 



t 



(N.2) 



where the function is tabulated in Table 4 for ma /nit < 4. For m^/m 4 > 4 we use the 
expansion Q 



t ' 144 



[311 + 24vr 2 + 282 In r + 90 In 2 r - 4r (40 + 6vr 2 + 15 In r + 18 In 2 r) 



+ r 2 (24209 - 6000vr 2 - 45420 In r - 18000 In 2 r)l, 

100 V >b 



(N.3) 



where r = t/h. For m t = 175 GeV and m# = 300 GeV 
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(2) 



1.245 . 
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Appendix O. 
The corrections. 



The corrections 62 Vi ~ aa s arise from gluon exchanges in quark electroweak loops |p4| 
(see also |82[). For two generations of light quarks (q = u, d, s, c) we have 

%V m (t, h) = 2- [^(^^)(c 2 - S 2 ) lnc 2 ] = (*=&^)(- .377) (0.1) 

5«V A (t, h) = h) = 2- [|(^^)(c 2 - , 2 + f /)] = (^^)(1.750) (0.2) 

%(*,/*) = (0.3) 

The result of calculations for the third generation is obtained in the form of fairly com- 
plicated functions of the t-quark mass: 

5iv m (t, h) = \(^){tM^) + (i - y* 2 )^) + (\ ~ I s2 ) ^ 

- 4(1 - S l) x tF.Cj) - 4^(0)} (0.4) 

siv A (t, h) = %v„{t, h) = 1 - -Mh+ 

6 71 4t 4 4t 

+ (1 - |0 2 [^i(^) - ^i'(^)] + (\~ \* + ^ 4 ) - 4^(0)} (0.5) 

5\V R {t, h) = ^(^Wi(^) " ^(1) - 4^(0) + \ hit}, (0.6) 

= l + ^a.(Mz)lnt ( °- 7) 

Note that are independent of m#. The functions Vi(r), Ai(r) and Fi(x) have rather 
complex form and were calculated in ||54|| . Their expansions for sufficiently small values of 
arguments are (we have added cubic terms to the expansions from 



where 



f,W=«(3)-^^ + ^ + ... (0.8) 

91 4Q fi89 3389 

AAr) = [-6C(3) - 3C(2) + — ] + r[4C(3) - — ] + r 2 — + r 3 — + ... (0.9) 

v ; 1 sv ; sv ; 4 J L sv ; 18 J 405 7 x 25 x 27 K ' 

3 1 93 1 95 

Fi(x) = [-|C(3) - -C(2) + ^] + x[C(3) - -C(2) + 

1 25 1 5 

+ x 2 [-C(2) + — 1 + x 3 [— C(2) + — 1 + ... (O.10) 

L 8 SV; 3x64 J [ 30 SV ; 72 J v ; 

where C(2) = vr 2 /6,C(3) = 1.2020569.... 
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Adding up the contributions (|Q.4|) -( pl6|) and using the expansions (|0.8|) -( p.lO|) , we 
obtain (up to terms ~ ^/t 3 ): 



8\V m {t, h) 



,a s {m t 



71 



)[-2.86t + 0.46 kit- 1.540 



$V A {t, h) = ^{t, h) 
,a s (mt 



,a s (m t 



71 



)[-2.86t + 0.493 



%V R (t,h) 



71 



)[-2.86t + 0.22 In t- 1.513 



0.68 


0.21 


_ & s {m t ) 


t 


t 2 J 


71 


0.19 


0.05 


_ a s {m t ) 


t 


t 2 J 


71 


0.42 


0.08 


_ a s (m t ) 


t 


t 2 J 


71 



-11.67) (O.ll) 

-10.10) (0.12) 
-11.88). (0.13) 



These formulas hold for m t > rriz- In the region m t < mz we either set 5 l 2 Vi = or make 
use of the massless limit 5\Vi = ^S^Vi. In any case this region gives negligible contribution 
to the global fit. 
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Appendix P. 
The S$Vi corrections. 

In the second order of weak interactions quadratic dependence on the mass of the higgs 



boson is given by expressions 83 



a ,mi, 0.747 

hV m = ^Hr x — r- = °- 0011 R1 

24tt m z <r 

rv m 2 1 1QQ 

^ = ^ = _(_| )x _ = 0.0057 (P.2) 
The numerical evaluations above were made with m# = 300 GeV. 
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